Award  Number:  W81XWH-04-1-0336 


AD 


TITLE:  Regulation  of  Leukocyte  Infiltration  into  Ovarian  Cancer  by  Tumor-Stroma 
Interactions,  a  Microarray  View  of  Cancer  Microenvironment 


PRINCIPAL  INVESTIGATOR:  Izhak  Haviv,  Ph.D. 

Mark  Smyth 

Carleen  Cullinane,  Ph.D. 
Michael  Kershaw,  Ph.D. 
Sarah  Russell,  Ph.D. 


CONTRACTING  ORGANIZATION:  University  of  Melbourne 

Victoria  3010  Australia 


REPORT  DATE:  March  2007 

TYPE  OF  REPORT:  Final 


PREPARED  FOR:  U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


DISTRIBUTION  STATEMENT:  Approved  for  Public  Release; 

Distribution  Unlimited 


The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision 
unless  so  designated  by  other  documentation. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

01-03-2007  Final  1  MAR  2004  -  28  FEB  2007 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 


Regulation  of  Leukocyte  Infiltration  into  Ovarian  Cancer  by  Tumor-Stroma 
Interactions,  a  Microarray  View  of  Cancer  Microenvironment 


5b.  GRANT  NUMBER 

W81XWH-04- 1-0336 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHOR(S) 

Izhak  Haviv,  Ph.D.  Mark  Smyth,  Carleen  Cullinane,  Ph.D.  Michael  Kershaw,  Ph.D. 
Sarah  Russell,  Ph.D. 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


E-Mail:  izhak.haviv@petermac.org 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

University  of  Melbourne 
Victoria  3010  Australia 


8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 
U.S.  Army  Medical  Research  and  Materiel  Command 
Fort  Detrick,  Maryland  21702-5012 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  Public  Release;  Distribution  Unlimited 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 


14.  ABSTRACT 

Despite  compelling  cell  biological  studies  and  histopathological  observations  incriminating  stromal  cells  in  tumorigenesis,  our  knowledge  of 
the  genes  that  mediate  changes  in  the  tumor  microenvironment  and  interactions  among  various  cell  types  in  epithelial  cancer  and  their  role 
in  tumorigenesis  is  limited.  Here,  we  describe  a  comprehensive  molecular  characterization  of  stromal-epithelial  cell  interactions,  using 
microarray  analysis  of  co-cultured  cell  pairs.  We  further  show  that  these  gene  expression  changes  indeed  are  common  feature  of  in  vivo 
profiles  of  invasive  epithelial  cancers.  Using  In  situ  hybridization  and  laser  capture  microdissection,  we  confirmed  previous  assignment  of 
which  cells  express  which  genes  in  response  to  co-culture.  We  provide  proof  of  principle  evidence  that  the  in  vitro  system  recapitulates  some 
aspects  of  cancer  cell  biology  and  may  contribute  new  markers  to  the  invasive  events.  We  directly  demonstrate  that  events  in  the  interaction 
of  epithelial  cells  and  their  neighboring  fibroblasts  play  initial  role  in  recruiting  inflammatory  cells  and  provoking  immune  surveillance.  We 
found  the  Epithelial-to-Mesenchymal  Transition  of  the  cancer  cells  (shRNA-mediated  E-Cadherin  knock  down)  attenuates  the  co-culture 
reaction,  suggesting  a  mechanism  for  cancer  cell  evasion  of  immune  surveillance. 


15.  SUBJECT  TERMS 

No  subject  terms  provided. 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

USAMRMC 

a.  REPORT 

U 

b.  ABSTRACT 

U 

c.  THIS  PAGE 

U 

uu 

59 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


Table  of  Contents 


Introduction . 4-5 

Body . 6-21 

Key  Research  Accomplishments . 21 

Reportable  Outcomes . 21 

Conclusions . 21 

References . 21-23 

Appendices . 24-59 


Annual  report  04  1  0336 


5/09/2007 


4  of  57 


OC030008 

Report  (according  to  Statement  of  Work)  of:  “the  Regulation  of  leukocyte  infiltration  into  ovarian 
cancer  by  tumor-stroma  interactions;  a  microarray  view  of  cancer  microenvironment”. 

Introduction: 

In  normal  tissues,  a  basement  membrane  separates  epithelia  and  mesenchymal  elements  (collectively 
termed  stroma),  which  support  the  survival  of  epithelial  cells  and  maintain  tissue  architecture  and  cellular 
boundaries  \  Invasion  of  stroma  by  transformed  epithelial  cells,  a  hallmark  of  cancer,  transgresses  the  normal 

••23 

boundary  that  separates  these  cell  lineages  and  results  in  abnormal  heterotypic  cell-cell  interactions  ’  .  The 
neighbors  of  the  epithelial  cancer  cells,  blood  vessels  4"6,  infiltrating  immune  cells  7‘ 8,  and  fibroblasts  4’ 9,  all 
exert  distinguishable  responses  to  this  process  (angiogenesis,  inflammation,  fibrosis).  These  responses  affect 
cancer  progression  in  a  several  ways.  Some  critical  invasion-promoting  enzymes  that  remodel  the  extracellular 
matrix  (ECM)  are  secreted  by  stromal  elements  10~12.  The  fibroblasts  that  reside  in  the  cancer  microenvironment 
are  phenotypically  different  than  normal  fibroblasts,  as  has  been  shown  for  breast  and  prostate  cancer.  They 
lack  the  ability  to  attenuate  the  growth  of  neighboring  transformed  epithelial  cells  13,  and  have  gained  the 
ability  to  accelerate  cancer  progression  14’  15.  The  mechanisms  responsible  for  this  phenotypic  change  may 
include  an  inherent  genomic  predisposition  ’  ,  somatic  genomic  aberrations  in  the  fibroblasts  '  ,  telomere 
attrition  ,  methylation  status  of  some  genes  “  ,  extracellular  matrix  density  and  morphological  changes  such 
as  desmoplasia  24.  While  the  incidence  of  visible  changes  in  the  stroma,  such  as  desmoplasia  are  well  recorded 
24,  the  incidence  of  the  above  described  phenotypic  and  genotypic  changes  in  fibroblasts  and  how  they 
contribute  to  carcinoma  pathology  are  harder  to  assess  and  discern  from  one  another.  Concurrently  definitive 
changes  in  the  gene  expression  patterns  of  the  carcinoma  fibroblasts  have  been  observed,  but  the  underlying 
mechanism  is  unknown  . 

DNA  microarrays  are  widely  used  to  study  changes  in  gene  expression  in  tumors.  Whilst  it  is  clear  that 
different  cell  populations  differ  in  their  gene  expression  ,  it  is  well  appreciated  that  the  changes  observed  in 
RNA  isolated  from  whole  tumors  represents  the  summation  of  expression  of  all  cellular  subtypes  2S. 
Hypothetically,  genomic  predisposition,  somatic  mutations,  telomere  attrition,  and  methylations  changes  in  the 
fibroblast  should  result  in  unique  gene  expression  signatures.  Therefore,  the  presence  of  these  unique  signatures 
in  tumour  expression  profiles  ’  ’  ’  could  be  used  to  infer  the  underlying  fibroblast  biology.  The  ability  of 
the  cancer  fibroblast  to  promote  or  to  attenuate  cancer  growth  should  also  be  associated  with  a  distinct  gene 
expression  signature.  If  the  genes  unique  to  each  of  these  changes  were  identified,  they  could  allow  us  to  assess 
the  frequency  of  the  associated  biology  in  human  cancers. 

This  DoD  idea  development  project  involves  expression  profiling  of  interactions  between  ovarian  cancer 
epithelial  cells  and  fibroblasts  with  the  specific  aim  of  identifying  the  underlying  mechanisms  of  tumor-stroma 
interaction.  Through  comparing  the  expression  response  to  co-culture,  to  those  observed  in  hundreds  of  human 
cancer  specimens,  we  validate  the  physiological  relevance  of  our  observations  to  human  ovarian  cancers.  We 
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defined  an  epithelial-fibroblast  interaction  signature,  supporting  previous  observations  29-3 ' ,  and  show  that  this 
signature  is  the  product  of  reciprocal  exchange  of  soluble  factors.  Fibroblast-produced  IL8  and  EGF  induce 
epithelial  cells  to  produce  TNFa,  inducing  the  fibroblasts  to  mount  an  inflammatory  reaction,  through 
chemokines,  tissue  remodeling,  and  the  extracellular  matrix.  This  interaction  may  prove  to  be  the  primary 
trigger  of  the  innate  immune  response  (parallel  to  TLR-mediated  responses)  and  the  underlying  mechanism  of 
the  chronic  wound  healing  responses  in  tumor  tissues  .  In  this  report  we  show  in  situ  hybridization  stains  of 
tissue  sections  with  cancer  specimens  that  demonstrate  the  same  genes  that  fibroblasts  express  in  response  to 
co-culture  in  vitro,  are  expressed  by  tumor-juxtaposed  fibroblasts  in  vivo.  To  address  the  physiological  effect  of 
infiltrating  immune  cells  on  ovarian  cancer  progression  we  introduced  human  ovarian  cancer  infiltrating 
immune  cells  into  nude  mice  in  xenograft  assays  of  ovarian  cancer  growth.  So  far,  these  experiments  have 
failed,  due  to  graft  versus  host  reaction.  We  are  in  the  process  of  depleting  mouse  antigen  reactive-cells  from 
the  tumor  infiltrating  immune  cells. 
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Task  1.  To  Extend  an  analysis  of  co-culture  regulated  genes  in  other  appropriate  cell  combinations, 
including  cells  derived  from  different  ovarian  tumor  stroma  and  epithelium  (Months  1-6): 


a.  Prepare  cells  from  different  tumors  with  defined  histology  and  use  MACS  separator  and  EpCAM- 
linked  magnetic  beads  to  separate  the  cells.  Success  is  assessed  on  flow  cytometry  using  2G3, 
Laminin,  FAPa  and  CK7  markers.  On  the  array  level  extend  the  arrays  from  10.5  k  cDNA  arrays  to 
19k  oligoarrays  (Compugen). 

We  have  so  far  isolated  five  cell  cultures  separated  from  ovarian  cancers,  that  are  >90%  enriched  for  CK7 
in  the  epithelial  component,  and  >80%  FAPa  for  the  non-epithelial  component.  At  this  moment,  however,  we 
have  not  got  enough  of  these  cells  to  perform  the  array  experiment.  We  are  currently  immortalizing  (papilloma 
E6/E7)  a  fraction  of  these  cells,  so  we  can  compare  the  immortalized  and  primary  cells  respond  the  same  way 
to  co-culture,  so  we  can  standardize  an  optimized  method.  In  addition,  we  have  characterized  the  response  to 
co-culture  on  Affymetrix  U133+2  full  genome  arrays,  for  three  of  these  culture  preps  (Figure  1). 

It  is  still  possible  that  signaling  pathways  other  than  TNFa,  TGFp,  IL8  and  EGF  that  we  found  so  far  do 

operate  in  these  patients, 
which  may  agree  with 
previous  co-cultured  pairs 
we  found,  where  the 
response  we  characterized 
so  far  was  missing.  Another 
source  for  novel  signaling 
pathways  between  the 
cancer  cell  and  the  stroma, 
we  will  obtain  viable  cells 
from  the  peritoneum  and 
omentum,  to  test  whether 
the  metastatic  environment 
involves  a  different 
response.  In  order  to 
specifically  look  for  more 
interaction  signaling 

pathways  than  the,  we  will 
continue  to  collect  viable 
cells  from  ovarian  cancer 
tissue,  but  will  only 
perform  the  co-culture 
experiment  if  the  frozen 
tissue  from  this  patient 
lacks  the  expression  profile 
we  found  is  correlated  with 
our  signaling  in  the  in  vivo 
expression  profiles.  The 
assumption  is  that  the 
absence  of  the  correlated  genes  indicates  the  absence  of  this  interaction  of  the  tumor  cell  with  the  stroma 
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b.  Express  ectopic  CDH1  in  OvCa  lines  with  mesenchymal  morphology  (Hey  and  A2780),  or 
pSUPER-based  siRNA  knock  down  vectors  for  this  gene  in  epithelial-like  OvCa  lines  (such  as 
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OVCAR-3).  Use  parental  and  derived  cells  in  co-culture  experiments.  Score  for  gene  expression 
changes  by  RT-PCR  (Months  2-6). 

A  common  reaction  to  the  co-culture  (like  2008),  involved  induced  expression  of  SNAI1  and  SNAI2,  or 
TWIST  1.  This  implies  that  the  interaction  with  the  theca  cells  induces  EMT  in  epithelial-like  ovarian  cancer 
cell  lines.  Similar  results  were  observed  with  the  fresh  cultures.  Activin  A  receptor  2,  JAK1  and  slit2  were 
elevated  in  the  non-reactive  lines,  while  STAT3  was  elevated  in  the  reactive  lines.  These  genes  were  still 
suggestive  of  cell  adhesion  and  polarity  as  playing  a  role,  in  a  similar  manner  as  proposed  in  the  application, 
but  not  necessarily  through  CDH1.  Furthermore,  since  2008  cells,  which  express  similar  amounts  of  CDH1  to 
OVCAR-3,  responds  to  co-culture  acutely,  had  segregated  between  the  minimal  response  to  co-culture  of 
OVCAR-3  cell  line,  and  its  strong  epithelial  morphology.  We  have  therefore  introduced  full  length  E-Cadherin 
plasmid  (Clone  ID  3546855,  OpenBiosystems®),  and  shRNA  plasmid,  targeting  e-Cadherin  (Clone  ID 
RHS41 86-97324246,  OpenBiosystems®),  to  OVCAR-3,  2008,  Hey  and  A2780  ovarian  cancer  cell  lines,  and 
subjected  the  modified  cells  to  a  co-culture  reaction.  This  part  of  the  work  is  not  concluded,  as  some  of  the  co¬ 
culture  combinations  need  to  be  repeated  for  result  confirmation.  However,  the  promising  and  surprising  result 
that  was  confirmed,  is  that  down-regulation  of  e-Cadherin  expression  resulted  in  robust  attenuation  of  co¬ 
culture  reaction  with  2008  cells.  While  this  results  runs  contrary  to  our  SOW  and  initial  application,  it  was 
reproduced,  and  we  are  certain  that  these  results  are  true.  This  means  that  epithelial  to  mesenchymal  transition, 
while  boosting  migratory  features  in  the  epithelial  cells,  also  turns  these  cells  into  stealth  cells,  in  terms  of  the 
co-culture  reaction.  The  consequence  of  this  change  in  2008  as  opposed  to  OVCAR-3  (control),  in  terms  of 
growth  in  xenograft  with  human  fibroblast  cell,  is  being  examined  currently. 

c.  Similar  to  Task  lb,  explore  defined  cellular  changes  to  the  ovarian  epithelial  cell  line  including 
DNA  damage  and  anoikis,  hormone  responses,  etc.  (Months  6-9). 

DNA  damage:  We  have  profiled  the  response  of  three  epithelial  cell  lines;  2780,  2008  and  OVCAR-3  to 
Cisplatinum  and  UV  treatments,  when  these  grow  alone,  or  in  a  tissue  culture  insert,  with  fibroblast  on  the 
bottom  chamber.  The  following  co-culture  responsive  genes  were  synergistically  induced  by  the  double 
treatment  of  co=culture  plus  DNA  damage;  MMP1,  MMP3,  MMP10,  FGF2,  COL5A2.  This  finding  is 
interesting  in  the  context  of  Col6A3  conferring  cisplatinum  resistance,  suggesting  that  cases  with  robust 
fibroblast  infiltrates  will  exhibit  worse  clinical  response  to  similar  treatment  (See  results  of  AOCS).  By 
contrast,  anoikis  and  hormone  treatment  did  not  change  the  course  of  the  co-culture  response  significantly. 

d.  Identify  potential  new  candidate  soluble  molecules  as  mediators  of  the  stroma  inflammation,  in  case 
the  expression  of  their  genes  are  affected  by  the  EMT.  Explore  those  factors  as  in  Task  2. 

The  factors  that  control  E-Cadherin  downregulation  are  still  being  explored,  however,  conditioned  media  from 
carcinoma  associated  fibroblasts  secrete  200  fold  more  IL6  than  normal  fibroblasts  in  co-culture.  In  our  hands, 
co-culture  in  the  presence  of  IL6  is  largely  attenuated,  supporting  the  possibility  that  EMT  quenches  co-culture 
response  and  implicating  IL6  in  common  EMT  in  breast  and  ovarian  cancer. 

e.  Microarray  profile  the  gene  expression  of  cultures  of  non-fibroblast  stromal  cells,  when  grown  in 
serum  free  medium  versus  conditioned  media  conditioned  by  epithelial  or  fibroblast  cultures,  as  well 
as  co-cultures  thereof.  Tested  cells  will  be  endothelial  (HUVEC),  Tumor  Infdtrated  Lymphocytes 
(TIL),  macrophages,  Theca  albea,  and  Granulosa  cells  (Months  6-12).  In  case  the  soluble  molecules 
appear  labile,  we  will  use  tissue  culture  inserts  (as  in  figure  5)  rather  than  conditioned  media. 

The  following  genes  were  expressed  at  differential  levels  upon  exposure  of  CD45+  cells  from  ovarian  cancer 
tissue,  in  response  to  epithelial-fibroblast  co-culture  medium:  CYBB,  MMP9,  IL4,  IL2RG,  CCL3,  PRG2,  LIF, 
CXCL9.  The  effects  of  these  changes  are  being  explored  in  the  proposed  xenograft  models  (not  completed).  We 
did 

f.  Define  binary  and  ternary  cell-cell  interactions  using  the  Task  Id.  Construct  a  table  of  inductive 
gene  expression  related  interactions  among  the  different  stroma  components. 

This  task  is  being  progressed  in  an  ongoing  fashion.  We  have  used  such  a  diagram  in  the  preparation  of  a 
successful  OC050265  grant  this  year.  This  diagram  will  continue  to  change  as  more  results  are  generated.  The 
hypotheses  that  such  a  diagram  is  generating  are  constantly  checked  against  expression  profiles  from  a 
constantly  growing  dataset  of  ovarian  cancer  expression  profiles,  as  well  as  xenograft  models  of  ovarian  cancer 
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tumour  growth.  We  include  a  number  of  such  co-transplanted  ovarian  cancer  xenograft  assays.  These  show  that 
the  effect  of  the  fibroblasts  is  constantly  increasing  the  cancer  growth.  All  the  cell  pairs  that  increase  in  growth 
in  response  to  co-transplantation,  share  a  gene  expression  response  to  co-culture,  which  is  similar  to  the  fresh 
co-culture  response,  as  shown  in  figure  1.  Also  integrated  are  expression  profiles  of  the  AOCS  ovarian  cancer 
patient  specimens.  These  provide  evidence  for  correlation  between  any  response  of  the  epithelial-fibroblast  co¬ 
cultures  and  the  infiltration  in  vivo  of  additional  cell  types  listed  in  Task  le  (using  their  cell  marker  mRNA  and 
expression  correlation  tools). 


Task  2.  To  identify  the  molecule(s)  that  mediate  the  co-culture  response. 

a.  Assign  responder-secretor  role  to  the  cells  in  the  co-culture.  For  example  in  point  a.2.7  epithelial 
cells  is  the  secretor  of  TNFa.  The  first  interaction  we  will  focus  on  is  the  response  of  the 
macrophages  to  a  fibroblast-derived  factor,  that  potentiates  the  macrophage  to  recruit  T  cells.  Here 
the  responder  is  the  macrophage  and  the  secretor  is  the  inflamed  fibroblast. 

In  task  1  we  identified  responses  of  leukocytes  to  the  co-culture,  and  identified  the  factors  induced  by  the 
binary  co-culture  with  the  fibroblasts,  mainly  chemokines,  that  would  induce  a  similar  response  in  the 
macrophage  and  natural  killer  cells,  as  the  co-culture  conditioned  media  would.  Interestingly,  the  chemokines 
and  their  corresponding  response  genes  in  the  leukocytes  are  correlated  in  ovarian  cancer  expression  profiles. 
We  are  now  characterizing  antibodies,  antagonists,  and  mainly  interfering  RNA  constructs,  in  attempt  to 
abrogate  the  leukocyte  response  to  co-culture.  This  part  of  the  work  is  done  in  collaboration  with  Shaun  McCall 
in  Adelaide  University,  a  world  expert  in  chemokines  signaling,  and  a  lab  loaded  with  chemokines  reagents. 
The  combinations  tested  are  mainly  focused  on  chemokines/cell  pairs  that  lead  to  migrational  effects.  We  are 
currently  troubleshooting  the  difficulty  to  attain  good  quality  array  results  from  the  leukocyte  population  (with 
and  without  co-culture  or  in  response  to  commercial  chemokines/cytokine),  as  our  responding  genes  were 
mainly  observed  with  real  time  PCR  on  candidates,  chosen  from  the  ovarian  cancer  expression  profiles  and  the 
literature,  such  as  CCL18  for  macrophages. 

b.  Construct  siRNA  vectors  (pRETROSUPER)  for  knockdown  of  candidate  molecules  (if  ovarian 
epithelial  cell  is  the  “responder”  in  the  co-culture,  siRNA  of  the  receptor  for  candidate  molecule 
will  be  targeted)  and  transfect  them  into  the  ovarian  epithelial  cells,  clone  plasmid-carrying 
derivative  cells.  Validate  that  these  vectors  mediate  expression  change  of  their  specific  genes  in 
stable  ovarian  epithelial  cell  clones  (Months  7-10). 

We  are  in  the  process  of  obtaining  and  managing  a  clone  set  (similar  to  cDNA  arrays)  of  the  Elledge/Hannon 
genome-wide  shRNA  lentivirus  library33.  This  will  spare  us  the  need  to  clone  and  construct  the  reagents  against 
the  target  gene  of  interest.  We  so  far  have  obtained  three  siRNA  vectors  that  modify  the  co-culture  reaction  in  a 
mnner  consistent  with  the  equivalent  siRNA  molecules  from  Ambion®.  Cells  modified  by  these  constructs  are 
currently  introduced  into  mice  xenografts,  with  the  aim  of  identifying  the  part  of  the  co-culture  response  that 
affects  cancer  growth,  and  thus  deserves  further  targeting. 

c.  Repeat  this  analysis  for  epithelial  response  to  fibroblast  that  appears  like  an  interferon-y  response 
(fibroblast  is  secretor  and  ovarian  epithelial  cell  is  responder),  as  in  Task  2a-e  (second  year). 

Neither  interferon  kappa,  or  combinations  of  the  factors  we  have  identified  to  accumulate  in  the  co-culture 
media,  such  as  TNFa,  IL8  and  IL6,  had  successfully  reproduced  the  interferon  response  we  observed  in  the  co¬ 
culture.  Since  JAK1  kinase  inhibitor,  and  siRNA  knocking  down  of  uPAR,  and  transfection  with  L1_VRAS,  all 
lead  to  complete  attenuation  of  the  signal  in  the  epithelial  cell,  we  suspect  that  the  interferon  response  is  a  late 
effect  of  the  RAS  signaling  induced  by  the  EGF-like  ligands  as  suggested  in  our  initial  analysis.  We  are 
exploring  this  possibility  now. 
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d.  If  Task  la-e  identify  new  co-culture  responses  that  are  robust  and  hint  at  clinically  relevant 
processes  in  the  AOCS  data,  repeat  this  analysis  for  this  new  response,  as  in  Task  2a-e  (year  3). 

Analysis  of  the  time  course  of  co-culture  found  TNFa  is  induced  in  8  hours  of  the  co-culture.  We  therefore 
tested  which  cytokines  induce  TNFa  expression  in  epithelial  cells.  We  found  EGF+IL8  or  EGF+LIF  induce  a 
200  fold  increase  in  TNFa  in  the  epithelial  cells.  Indeed,  siRNA  against  IL8RA  attenuates  the  co-culture  (not 
shown).  This  observation  is  consistent  with  IL8  being  expressed  by  the  fibroblasts.  We  are  testing  siRNA 
modified  fibroblasts  and  epithelial  cells  in  ovarian  cancer  xenograft,  to  see  if  (1)  the  co-culture  signature  is 
modified,  or  (2)  cancer  growth  is  modified,  if  TNFa  signaling  is  impeded  in  the  fibroblast,  or  if  IL8  signaling  is 
impeded  in  the  epithelial  cells. 

e.  Acquire  neutralizing  antibodies  against  candidate  molecule  (that  fulfilled  Task  2b  requirements). 
Apply  these  antibodies  onto  a  co-culture  experiment.  Check  if  antibody  diminishes  the  co-culture 
response  (Month  10). 

f.  Construct  siRNA  vectors  (pRETROSUPER)  for  knockdown  of  candidate  molecules  (if  ovarian 
epithelial  cell  is  the  “responder”  in  the  co-culture,  siRNA  of  the  receptor  for  candidate  molecule 
will  be  targeted)  and  transfect  them  into  the  ovarian  epithelial  cells,  clone  plasmid-carrying 
derivative  cells.  Validate  that  these  vectors  mediate  expression  change  of  their  specific  genes  in 
stable  ovarian  epithelial  cell  clones  (Months  7-10). 

g.  Confirm  the  role  of  the  candidate  molecule  in  the  co-culture.  Perform  co-culture  microarrays  on 
derived  clones,  compare  to  parental  lines  for  co-culture  response  (Months  8-12). 

We  are  in  the  process  of  obtaining  and  managing  a  clone  set  (similar  to  cDNA  arrays)  of  the  Elledge/Hannon 
genome-wide  shRNA  lentivirus  library33.  This  will  spare  us  the  need  to  clone  and  construct  the  reagents  against 
the  target  gene  of  interest.  Clones  carrying  these  clones  were  derived  from  2008,  OVCAR-3,  and  fibroblasts, 
and  the  effect  of  each  knock-down  on  the  co-culture  was  tested  by  a  95  real  time  PCR  assays  from  cDNA 
derived  of  epithelial-fibroblast  co-culture  experiments  using  the  modified  cells  with  normal  counterpart.  We  so 
far  have  obtained  three  siRNA  vectors  that  modify  the  co-culture  reaction  in  a  manner  consistent  with  the 
equivalent  siRNA  molecules  from  Ambion®.  Cells  modified  by  these  constructs  are  currently  introduced  into 
mice  xenografts,  with  the  aim  of  identifying  the  part  of  the  co-culture  response  that  affects  cancer  growth,  and 
thus  deserves  further  targeting. 

A  partial  list  of  the  tested  shRNA  constructs  in  co-culture#: 


# 

SYMBOL 

CLONE  ID 

Library  clone 

\co-culture  Effect* 

1 

AEBP1 

RHS3979-9588632 

TRCN0000021 224 

3F 

2 

CDH1 

RHS41 86-97324246 

VLHSJ00 

5E 

3 

CDH11 

RHS3979-9621532 

TRCN0000054333 

1 _ | 

4 

COLEC11 

RHS1 764-9680909 

V2HS_98658 

1 _ |  1 

5 

CTSK 

RHS41 86-97557536 

VLHS_77990 

1 

6 

DDR2 

RHS3979-9569897 

TRCN0000001417 

3E 

7 

FAP 

RHS41 86-975571 99 

VLHS_74865 

1 

8 

FN1 

RHS41 86-973241 96 

VLHS_99078 

1 

9 

HGF 

RHS41 86-97322497 

VLHS_41775 

1 

10 

IL8 

RHS41 86-97323076 

VLHS_96688 

4E 

11 

INHBA 

RHS3979-9626447 

TRCN0000059263 

1 

12 

IRF7 

RHS3979-9582548 

TRCN00000 14859 

2E 

13 

PDGFD 

RHS41 86-97321 033 

VLHS_1 22453 

1 _ | 

14 

PDGFRB 

RHS41 86-97551 128 

VLHS_154717 

4F  1 

15 

PLAUR 

RHS41 86-9755451 5 

VLHS_1 55000 

4E 
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16 

RAB31  RHS3979-96 14932  TRCN0000047733  1 

17 

RON  RHS3979-9571732  TRCN0000003106  1 

18 

RUNX1  RHS41 86-97658054  VLHS  135256  1 

19 

SFRP4  RHS3979-95821 77  TRCN0000014488  1 

20 

SMAD4  J  RHS41 86-97548297  VLHS  22463  |1_ 

21 

SMAD7  J  RHS3979-9586792  TRCN0000019384  1 

22 

Snail  RHS41 86-9754931 4  VLHS_1 83302  3E 

23 

SNAI2  RHS41 86-97323250  VLHS_1 38122  3E 

24 

SPARC  RHS3979-9578042  TRCN0000008709  1 

25 

TGFBRII  RHS3979-97079721  TRCN00000 10444  1 

26 

Tnfrl  RHS41 86-97550627  VLHS_79261  5F/3E 

27 

CCL5  RHS1 764-9688686  V2HSJ31260  1 

28 

CCL5  RHS1 764-9687862  V2HS  31263  1 

29 

CCL7  J  RHS41 86-97323992  VLHS  42231  |1_ 

30 

CCL8  J  RHS41 86-97551 308  VLHS  42021  |l_ 

31 

CCL1 1  J  RHS1 764-9399657  V2HS  31271  |l_ 

32 

CCL11  RHS1 764-9209707  V2HS_31275  1 

33 

CCL11  RHS1 764-9690623  V2HS_31274  1 

34 

CXCL6  RHS41 86-976571 79  VLHS_1 6254  1 

35 

CXCL12  RHS41 86-97555397  VLHS_61280  1 

36 

CXCL12  J  RHS41 86-97554397  VLHS  61277  |1_ 

37 

CXCL12  RHS41 86-97322255  VLHS_96819  1 

38 

CXCL12  RHS41 86-97553733  VLHS  80578  1 

39 

CXCL12  RHS41 86-975581 30  VLHS  214349  T~ 

40 

ITGB1  RHS41 86-97550631  VLHS_1 18535  1 

41 

ITGB1  J  RHS41 86-97323883  VLHS  219947  |1_ 

42 

ITGB1  RHS41 86-97549749  VLHS_1 18533  1 

43 

NFKBIA  RHS41 86-97552624  VLHS_51129  4E 

44 

NFKBIE  RHS3979-962401 7  TRCN0000056833  4E 

45 

NFKBIE  RHS3979-962401 8  TRCN0000056834  3E 

46 

NFKBIE  RHS3979-962401 9  TRCN0000056835  \2E 

47 

NFKBIE  RHS3979-9624021  TRCN0000056837  3E 

48 

LYN  RHS41 86-97658662"  VLHS  119208  |1_ 

49 

LYN  J  RHS41 86-97324707  VLHS  119204  |1_ 

50 

IFITM2  RHS41 86-97553774  VLHS_234126  1 

51 

NMI  RHS1 764-9499660  V2HS_41117  1 

52 

AIF1  RHS1 764-9681 233  V2HS_1 32490  1 

#-These  experiments  have  not  been  reproduced  more  than  once  so  far,  and  are  major  cause  of  the  delay  in 
reporting. 
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*-Results  presented  in  1-5,  where  1  is  no  effect  on  co-culture  response,  and  5  is  total  shut-down  of  the  response. 
E  following  the  number  signifies  that  the  knock  down  which  affected  co-culture  response  occurred  in  the 
epithelial  cells,  and  F  represents  knock  down  in  fibroblasts. 

h.  Repeat  this  analysis  for  epithelial  response  to  fibroblast  that  appears  like  an  interferon-y  response 
(fibroblast  is  secretor  and  ovarian  epithelial  cell  is  responder),  as  in  Task  2a-e  (second  year). 

i.  If  Task  la-e  identify  new  co-culture  responses  that  are  robust  and  hint  at  clinically  relevant 
processes  in  the  AOCS  data,  repeat  this  analysis  for  this  new  response,  as  in  Task  2a-e  (year  3). 

The  most  significant  observation  from  the  shRNA-mediated  knock  down  experiments  is  the  effect  of 
CDHl/E-cadherin.  The  expression  of  this  gene  in  AOCS  cases  is  robustly  different  among  the  molecular 
classes  of  OvCa,  with  the  highest  expression  in  class  5,  which  is  lacking  both  the  infiltrating  inflammatory 
signature,  as  well  as  the  fibroblast  signature.  Xenograft  model  experiment  with  CDH1  knock-down  are 
ongoing. 

Task  5.  To  compare  the  expression  pattern  of  key  molecules  that  come  up  in  co-culture  in  surgical 
specimens  from  CDRP  project  OC000109,  with  the  aim  of  assessing  how  general  our  observations  are, 
and  in  order  to  relate  our  observations  with  in  vivo  clinical  features  such  as  lymphocyte  infiltration  in 
tumors  and  patient  survival  (parallel  to  Task  2,  after  each  candidate  is  identified). 

a.  Establish  as  many  as  possible  Tissue  arrays  with  needle  punch  sections  arrayed  on  glass  slides. 

255  of  AOCS  ovarian  cancers  and  parallel  tissue  banks  have  been  generated  in  Peter  MacCallum  Cancer  Centre 
for  this  project. 

b.  Obtain  antibodies  against  candidate  molecules  that  are  reactive  and  specific  in 
Immunohistochemistry. 

A  drawback  in  the  task  is  that  many  of  the  interesting  genes  are  not  reactive  on  formalin  fixed  material  and 
require  fresh  frozen  tissue.  This  means  we  stain  individual  specimen,  rather  than  the  high  throughput  of  the 
tissue  array.  As  an  alternative  we  resort  to  Task  3c,  which  means  ISH  and  we  have  the  assay  pipelined  and 
validated  many  of  the  CRGs  in  vivo  (see  below). 

c.  Stain  for  candidate  molecules  by  both  in  situ  hybridization  for  RNA  expression,  as  well  as 
immunohistochemistry  to  assess  correlation  of  expression  with  biological  processes  that  are  the 
focus  of  this  study,  such  as  TIL  (CD3+  cells  in  the  tumor  sections),  angiogenesis,  inflammatory 
reactions,  invasion,  and  apoptosis. 

AOCS  tissue  arrays  were  stained  with  CRG-specific  antisense  probes  for  in  situ  hybridization.  Sense  probe 
showed  no  staining.  Shown  below  are  genes  that  our  previous  analysis  assigned  as  fibroblast  response  genes. 
Indeed  in  vivo  sections  show  that  the  genes  are  expressed  in  fibroblasts  immediately  juxtaposed  to  cancer  cells. 
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Figure  1 :  In  situ  hybridization  of  CRGs  on  cancer  specimen  sections,  shown  in  clockwise  sequence  from  top  left;  MMP3, 
TNFAIP6,  INHBA,  SPARC  and  PDGFRB.  Similar  results  were  obtained  from  TNC,  SERPINE1  (PAI1),  IL8  and  OAS1. 
In  all  cases,  the  staining  cells  are  almost  exclusively  the  fibroblasts,  immediately  juxtaposed  to  cancer  cells. 
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d.  Assess  our  hypotheses  in  light  of  the  expression  associations  with  the  known  clinical  data  of  those 
specimens. 


Immune  Cell  Signature  in  Absence  of  Myofibroblast 
Signature  Correlates  with  Better  Survival 

K3  K4  Progression 

Poor  Good  Free  Survival 


Overall 

Survival 


AOCS  ovarian  cancer  specimens  subdivide  into  six  molecular  classes,  the  two  most  significantly  distinct  in 
terms  of  clinical  data  share  infiltrating  inflammatory  signature,  and  differ  in  (activated/myo)  fibroblastic 
signature.  LCM  of  six  such  samples  demonstrated  that  all  genes  in  these  signatures  are  expressed  in  the  stromal 
infiltrating  cells. 
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Laser  Capture  Microdissection  Validates  Stroma  Specific  Gene  Expression 

G1 

Stroma 
Specific 
Genes 


5  Samples  (G1)  Tumour  and  Stroma 
NuGEN  Ovation  Pico  Amplification 
Affymetrix  U133  Plus  2  arrays 
Paired  sample  SAM  Tumour  vs  Stroma 
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Task  4.  To  characterize  the  biological  significance  of  the  gene  expression  changes  observed  in  co-cultures 
and  their  effects  on  cancer  progression  and  immune  surveillance  of  cancers  in  vitro  and  in  vivo. 

a.  Use  flow  cytometry  kits  like  APO-BrdU™  in  CytoPerm-Cytofix  ™  buffers,  for  the  analysis  of  cell 
cycle  and  cell  death  process  in  cells  growing  in  monocultures  versus  co-cultures.  Individual  cell 
populations  are  analyzed  in  the  mixed  cell  populations  by  gating  the  results  with  CK7  as  marker 
for  epithelial  cells. 

b.  Construct  Ectopic  expression  vectors  for  molecules  that  are  found  in  Task  2  to  confer  the  co¬ 
culture  responses.  Use  a  modified  version  of  pBABE  vectors,  carrying  ovarian  3  specific  promoter. 

c.  Construct  Ectopic  expression  vectors  for  molecules  that  are  found  in  Task  2  to  confer  the  co¬ 
culture  responses.  Use  a  modified  version  of  pBABE  vectors,  carrying  ovarian  3  specific  promoter. 

d.  Transfect  obtained  vector  into  epithelial  cancer  cell  lines  that  are  less  active  in  the  co-culture 
assays,  such  as  A2780  (Months  7-10). 

e.  Construct  knock  down  vectors  siRNA  vectors  (pRETROSUPER  vector),  for  the  same  genes  as  in 
task  4b. 

f.  Transfect  those  into  epithelial  cancer  cell  lines  that  are  less  active  in  the  co-culture  assays,  such  as 
2008  or  Hey. 

g.  Construct  expression  vectors  for  negative  dominant  genes  such  as  SOCS  to  disrupt  the  interferon 
response  in  the  epithelial  cells. 

h.  Transfect  epithelial  cancer  cell  line  A2780  and  2008  with  these  vectors. 

i.  Perform  co-culture  microarrays  on  derived  clones,  compare  to  parental  lines  for  co-culture 
response  (Months  8-12). 

j .  Choose  cell  derivatives  (carrying  a  gene  expression  aberration  for  the  molecules  mediating  the  co¬ 
culture  response),  and  characterize  their  malignant  potential  in  mouse  xenograft  models. 
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k.  Assess  the  response  to  the  transplant  by  the  mouse  stroma  (compare  gene-modified  cells  to 
parental). 

Experiment  la  (STR:01a) 

5  x  106  2008  cells  +/-  5  x  105  FF  cells 


Experiment  lb  (STR:01b) 

3  xlO6  A2780  cells  +  5  x  105  FF  cells 

No  growth  of  A2780  cells  alone;  graph  represents  growth  of  A2780  +  FF  cells 


A2780  cells  +  FF 
A2780 
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Experiment  4  (STR:04) 

3.5  x  106  CH-1  cells  +  5  x  105  FF  cells 

No  growth  of  CHI  cells  alone.  Only  one  tumour  in  the  CHI  +  FF  group  grew. 
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Experiment  5  (STR:05) 


3  x  106  SKOV-3  cells  +/-  5  x  104  FF  cells 


^  400 
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Experiment  6  (STR:06) 

5  x  106  DOV-13  cells  +/-  5  x  105  FF  cells 


0  5  10  15  20  25  30  35  40  45 

Days 


Experiment  9  (STR:09) 

3  x  106  SKOV-3  cells  +/-  5  x  10s  FF  cells 

Tumours  shown  for  individual  mice-  some  tumours  from  each  group  did  not  take. 
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1.  Introduce  human  stromal  components  together  with  the  ovarian  cancer  cell  line.  Start  with 
fibroblasts,  then  macrophages. 

Open  Biosystem  ®  TRC  and  Elledge/Hannon  constructs  of  miRNA  in  lentiviruses,  directed  against  IL8RA, 
JNK,  STAT3,  IicKy  and  TNFR1  are  all  stably  transfected  into  fibroblasts  and  ovarian  cancer  cells.  Knocked 
down  colonies  are  being  introduced  into  xenograft  cancer  model,  to  test  the  degree  of  synergy  with  the 
fibroblasts,  in  conditions,  which  abrogate  the  interaction  we  characterized  in  vitro.  Growth  curves  will  be 
performed  as  described  in  the  previous  report. 


m.  For  Fymphocyte  recruitment,  ovarian  epithelial  cancer  cells  will  be  used  from  our  own  tissue  bank 
collection.  These  cells  have  corresponding  peripheral  blood  lymphocytes,  which  avoids  allograft 
rejection  reactions  in  the  xenograft.  Obtain  these  lymphocytes,  co-inject  them  with  the  cancer 
cells,  fibroblasts,  and  macrophages,  all  from  the  same  patient  (so  far  have  such  cell  cultures  from 
two  patients). 

Most  xenografts  performed  this  year  focused  on  mice  carrying  human  leukocytes.  We  found  that  in  the  absence 
of  immune  compromising  drugs,  the  mice  all  die  very  quickly  (within  ten  days).  We  are  currently  trying  to 
rectify  this  problem  by  depleting  mouse  reactive  cells  from  our  leukocyte  and  TIF  populations.  If  these  immune 
cells  become  stable  in  the  mice,  we  will  proceed  with  this  proposed  experiment.  Otherwise,  we  will  resort  to 
attenuating  the  immune  response  with  drugs  that  prolong  graft  survival,  which  would  be  unfortunate  to  the  aim 
of  the  experiment.  However,  in  vitro,  we  have  characterized  the  migration  of  CD45+  cells  to  conditioned 
medium  of  co-cultures  and  to  defined  cytokines  and  chemokines. 


Migration  of  PBLstoCCL7  cfto 
migration  to  control 


Time  points  at  which  PBLswere  taken  and 
counted 


n.  Follow  mice  for  survival,  tumor  infiltration  with  lymphocytes,  inflammatory  reaction,  cancer 
histology,  tumor  burden  time  course,  immunohistochemistry  with  antibodies  against  co-culture 
induced  gene  products  (such  as  CXCF9  and  MMP2). 

Experiments  are  ongoing. 
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Key  Research  Accomplishments: 

1)  Identification  of  genes  that  best  correlate  with  dynamic  interactions  between  epithelial  cancer  cells  and 
neighboring  fibroblasts. 

2)  Identification  of  the  cytokines  that  are  producing  this  gene  expression  response. 

3)  Identification  of  intrinsic  changes  in  the  cancer  cell  that  modulate  the  extent  of  interaction  with 
fibroblasts,  and  consequent  inflammatory  response. 

4)  Linking  these  interaction  with  neighboring  fibroblasts  with  ovarian  cancer  patient  outcome,  the  novelty 
of  this  finding  is  the  emphasis  on  the  state  of  activation  of  the  fibroblasts,  as  opposed  to  previously 
reported  relative  infiltration  and  positioning  of  T-cell  populations. 

Reportable  Outcomes: 


Enclosed  please  find  a  final  version  of  paper,  in  final  stages  of  preparation  before  submission.  This  is  where 
most  of  the  effort  is  going  currently,  as  this  work  is  long  due,  and  DoD  deserves  acknowledgement  for  its 
support. 

Conclusions: 

The  results  presented  in  this  report  suggest  that  inhibiting  the  activation  of  fibroblasts  in  tumor 
environment  may  have  clinical  benefit  to  ovarian  epithelial  cancer  patient  outcome.  This  may  be  done  by  anti¬ 
inflammatory  reagents,  such  as  TNF-neutralizing  antibodies,  as  we  used,  or  by  inhibitors  of  receptor  tyrosine 
kinases,  which  are  vital  to  the  activated  fibroblasts,  such  as  PDGFRB  and  FGFR2  (see  shRNA  results). 
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Gene  expression  analysis  of  interacting  fibroblast  and  epithelial  ovarian  tumor  cells 
recapitulates  a  cell  inflammatory  signature  common  in  cancer 
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Abstract:  Gene  expression  analysis  of  a  whole  tumor  sums  the  mRNA  contributions  of 
cell  lineage,  cell-autonomous  genomic  changes,  and  non-cell  autonomous  cell-matrix  and 
cell-cell  interactions.  Extrapolating  the  latter  from  the  reaction  of  cells  to  co-culture  was 
used  to  derive  intercellular  signals  that  contribute  to  cross-talk  between  ovarian  epithelial 
cells  and  fibroblasts.  A  reciprocal  interaction,  where  fibroblast-derived  IL8  and  EGF 
provoke  the  epithelial  cells  to  express  TNFa,  thus  driving  inflammatory  gene  products  in 
the  fibroblasts,  was  found  to  be  a  highly  consistent  feature  of  human  cancer  expression 
profiles.  This  work  offers  a  compendium  of  experimental  growth  conditions  and  their 
consequential  gene  expression  changes  which  may  be  suitable  for  future  attempts  to 
modulate  the  cancer  microenvironment,  by  focusing  on  the  signals  that  are  most  robust 
and  consistent  in  a  large  fraction  of  cancer  patients  of  a  variety  of  carcinomas. 

Running  Title:  Reciprocal  interactions  between  fibroblastic  and  epithelial  cells  drive 
cancer  specific  expression  profiles. 


'Trescowthick  Research  Laboratories  and  the  Ian  Potter  Centre  for  Cancer  Genomics  and  Predictive 
Medicine.  Peter  MacCallum  Cancer  Centre,  St.  Andrews  Place,  East  Melbourne,  VIC  Australia  3002. 
2University  of  Melbourne,  Department  of  Biochemistry  and  Molecular  Biology,  Parkville,  Melbourne. 

3  National  Infomiation  and  Computing  Technology  Australia  Institute,  Australian  National  University, 
Canberra 

4  University  of  Melbourne,  Department  of  Medicine,  Parkville,  Melbourne. 

3  New  South  Wales  Pathology  Laboratory, 

6  Cambridge,  UK 

7  Walter  and  Eliza  Hall  Institute,  Parkville,  Vic  3052,  Australia 

g 

Corresponding  author,  d.bowtell@petermac.org 


1 


NOTE:  total  article  <  55,000  characters 

Introduction 

In  normal  tissues,  a  basement  membrane  separates  epithelial  and  mesenchymal  elements 
(collectively  tenned  stroma),  which  support  the  survival  of  epithelial  cells  and  maintain  tissue 
architecture  and  cellular  boundaries  (Jacks  and  Weinberg,  2002).  Invasion  of  stroma  by 
transformed  epithelial  cells,  a  hallmark  of  cancer,  transgresses  the  normal  boundary  that 
separates  these  cell  lineages  and  results  in  abnormal  heterotypic  cell-cell  interactions  (Liotta  et 
al.,  1980;  Yang  et  ah,  2002).  The  neighbors  of  the  epithelial  cancer  cells,  blood  vessels  (Allinen 
et  al.,  2004;  Bamberger  and  Perrett,  2002;  Zurita  et  al.,  2003),  infiltrating  immune  cells  (Haskill 
et  al.,  1982;  Zhang  et  al.,  2003),  and  fibroblasts  (Allinen  et  al.,  2004;  Kohno  et  al.,  1982),  all 
exert  distinguishable  responses  to  this  process  (angiogenesis,  inflammation,  fibrosis).  These 
responses  affect  cancer  progression  in  several  ways,  importantly,  critical  invasion-promoting 
enzymes  that  remodel  the  extracellular  matrix  (ECM)  are  secreted  by  stromal  elements 
(Hiratsuka  et  al.,  2002;  Lynch  and  Matrisian,  2002;  Sugiyama  et  al.,  2001).  The  fibroblasts  that 
reside  in  the  cancer  microenvironment  are  phenotypically  different  than  normal  fibroblasts,  as 
has  been  shown  for  breast  and  prostate  cancer.  They  lack  the  ability  to  attenuate  the  growth  of 
neighboring  transfonned  epithelial  cells  (Sadlonova  et  al.,  2005),  and  have  gained  the  ability  to 
accelerate  cancer  progression  (Olumi  et  al.,  1999;  Orimo  et  al.,  2005).  Several  mechanisms  may 
be  responsible  for  this  phenotypic  change  in  tumor  related  fibroblasts,  including,  an  inherent 
genomic  predisposition  (Deshpande  et  al.,  2005;  Morgan  et  al.,  2005),  somatic  genomic 
aberrations  in  the  fibroblasts  (Kurose  et  al.,  2002;  Moinfar  et  al.,  2000;  Tuhkanen  et  al.,  2004), 
telomere  attrition  (Parrinello  et  al.,  2005),  changes  in  the  methylation  status  of  some  genes  (Hu  et 
al.,  2005),  changes  of  the  extracellular  matrix  density  (Paszek  et  al.,  2005)  and  morphological 
alterations  such  as  desmoplasia  (Madri  and  Carter,  1984).  While  the  incidence  of  visible  changes 
in  the  stroma,  such  as  desmoplasia,  are  well  recorded  (Madri  and  Carter,  1984),  the  incidence  of 
the  above  described  phenotypic  behavior  and  genotypic  changes  in  fibroblasts  and  their  exact 
contribution  to  carcinoma  pathology  are  harder  to  assess  and  discern  from  one  another. 
Concurrently,  definitive  changes  in  the  gene  expression  patterns  of  the  carcinoma  fibroblasts 
have  been  observed,  but  the  inciting  mechanism  remains  unknown  (Allinen  et  al.,  2004; 
Cmogorac-Jurcevic  et  al.,  2001;  Hu  et  al.,  2005;  Parker  et  al.,  2005). 
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DNA  microarrays  are  widely  used  to  study  changes  of  gene  expression  in  tumors.  Whilst 
it  is  clear  that  different  cell  populations  differ  in  their  gene  expression  (Perou  et  ah,  1999),  it  is 
also  well  appreciated  that  the  changes  observed  in  RNA  isolated  from  whole  tumors  represent 
the  summation  of  expression  of  all  cellular  subtypes  (Stuart  et  al.,  2004).  Hypothetically, 
genomic  predisposition,  somatic  mutations,  telomere  attrition,  and  methylation  changes  in  the 
fibroblast  should  each  result  in  unique  gene  expression  signatures.  Therefore,  the  presence  of 
these  unique  signatures  in  tumor  expression  profiles  (Allinen  et  al.,  2004;  Cmogorac-Jurcevic  et 
al.,  2001;  Hu  et  al.,  2005;  Parker  et  al.,  2005)  could  be  used  to  infer  the  underlying  fibroblast 
biology.  Therefore,  the  ability  of  the  cancer  fibroblast  to  promote  or  to  attenuate  cancer  growth 
should  also  be  associated  with  a  distinct  gene  expression  signature.  If  the  genes  unique  to  each  of 
these  changes  were  identified,  they  could  allow  us  to  assess  the  frequency  of  the  associated 
fibroblast  biology  in  human  cancers. 

We  have  expression  profiled  tumor-stromal  interactions  between  cancer  epithelial  cells 
and  fibroblasts  with  the  specific  aim  of  identifying  the  underlying  mechanisms  of  tumor-stroma 
interaction,  and  have  compared  the  profiles  to  those  observed  in  hundreds  of  human  cancer 
specimens.  We  defined  an  epithelial-fibroblast  interaction  signature,  supporting  previous 
observations  (Fromigue  et  al.,  2003;  Gallagher  et  al.,  2005;  Sato  et  al.,  2004),  and  show  that  this 
signature  is  the  product  of  reciprocal  exchange  of  soluble  factors.  Fibroblast-produced  IL8  and 
EGF  induce  epithelial  cells  to  produce  TNFa,  inducing  the  fibroblasts  to  mount  an  inflammatory 
reaction,  through  chemokines,  tissue  remodeling,  and  the  extracellular  matrix.  This  interaction 
may  prove  to  be  the  primary  trigger  of  the  innate  immune  response  (parallel  to  toll  like  receptor- 
mediated  responses)  and  the  underlying  mechanism  of  the  chronic  wound  healing  responses  in 
tumor  tissues  (Dvorak,  1986). 

Results: 

Co-culture  of  ovarian  epithelial  tumor  cells  and  fibroblasts  results  in  the  induction  of 
specific  gene  expression 

To  better  understand  gene  expression  changes  that  are  attributable  to  intercellular  cross¬ 
talk,  we  compared  mono-cultured  cells,  with  co-cultured  cells,  by  microarray  gene  expression 
profiling  (Figure  1  A).  Equal  numbers  of  cells  were  grown  either  in  co-culture  or  in  mono-culture 
for  16-24  hours.  After  pooling  of  the  monoculture  cells,  the  RNA  was  isolated  from  both  the 
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monoculture  and  co-culture  cell  mixtures,  labeled  with  fluorescent  nucleotides  and  the  targets 
competitively  hybridized  with  cDNA  microarrays  (Fig.  1A).  If  there  was  no  intercellular 
communication  in  co-culture,  we  would  expect  no  change  in  gene  expression  (i.e.  no  log  ratios 
substantially  different  from  zero),  however,  we  found  that  the  expression  of  a  large  number  of 
genes  were  reproducibly  altered  by  co-culture  (Table  1,  column  M),  particularly  genes  involved 
in  ECM  production  and  remodeling,  cell  signaling  and  inflammation.  The  ratio  of  cell  types  in 
the  mixed  mono-cultures  and  co-cultures  were  comparable  based  on  FACS  and  cell-type 
expression  profiles  (supplement  1,  Supplementary  Figure  1A  &  B).  Figure  IB  shows  how  genes 
were  selected  for  further  analysis,  based  on  their  fold  expression  changes  and  the  consistency  of 
these  changes  across  experimental  repeats. 

Assessing  the  conservation  of  co-culture  induced  genes  across  multiple  cell  pairs. 

Several  co-culture  experiments  characterizing  the  fibroblast-epithelial  interaction  by  gene 
expression  analysis  have  been  perfonned  (Fromigue  et  al.,  2003;  Gallagher  et  al.,  2005).  In  each 
case  one  fibro-epithelial  cell  pair  was  used,  assuming  that  the  gene  expression  response  is 
generalizable.  We  have  used  multiple  fibro-epithelial  cell  pairs  to  assess  this  assumption.  The  co¬ 
culture  experiment  described  above  was  repeated  using  15  ovarian  cancer  cell  lines  and  one 
fibroblast  primary  cell  culture.  For  each  cell  pair,  co-culture  and  mixed  mono-cultures  cDNAs 
were  prepared  and  compared  as  in  Figure  1A.  We  performed  real  time  PCR  using  the  top  125 
genes  that  were  upregulated  in  the  original  experiment,  described  above  (OAS2  and  OAS3 
omitted,  see  Supplementary  Table  I).  Expression  level  ratios  of  each  gene  between  co-culture 
and  mixed-monocultures  were  normalized  (Vandesompele  et  al.,  2002)  and  used  to  produce  a 
principle  component  analysis,  measuring  the  degree  of  similarity  between  the  cell  pairs  (Figure 
2A),  and  to  generate  a  hierarchical  cluster  (Figure  2B).  The  overall  change  in  gene  expression  is 
also  represented  by  box  plots  (Figure  2C).  The  median  expression  change  in  co-culture  was 
higher  than  two  fold  for  12  of  the  cell  pairs  (Figure  2C  and  Supplementary  Table  II).  By  contrast, 
some  cell  pairs  exhibited  little  gene  expression  change  (Figure  2A  and  B,  blue).  The  56  most 
robust  and  consistently  upregulated  co-culture  genes  are  marked  in  bold  in  Table  1.  The  same 
genes  were  found  to  respond  to  co-culture  when  the  experiments  were  repeated  with  fibroblasts 
and  epithelial  cells  that  derived  from  disparate  organs,  such  as  skin,  breast,  ovarian,  gastric  and 
lung  (data  not  shown).  Of  particular  interest  is  the  observation  that  a  robust  response  also 
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occurred  when  non-transformed  epithelial  cells  were  co-cultured  with  fibroblasts  (e.g.  Moody 
human  ovarian  surface  epithelium,  Figure  2B,  and  other,  data  not  shown).  This  suggested  that  the 
phenomenon  we  are  modeling  in  the  in  vitro  co-culture  is  relatively  conserved  across  different 
types  of  carcinomas,  and  would  justify  analysis  of  these  genes  in  expression  profiles  of  multiple 
types  of  human  cancer  specimens. 

Co-culture  induces  genes  involved  in  tissue  remodeling,  inflammation,  and  interferon 
responses. 

One  hundred  and  twenty  seven  genes  were  identified  that  were  consistently  up  regulated 
by  co-culture  (Table  1,  column  M).  The  upregulated  genes  are  strikingly  enriched  for  molecules 
that  mediate  inflammatory  responses  and  tissue  remodeling,  including  13  inflammatory  mediator 
genes,  or  genes  reported  as  TNFa-responsive,  such  as  IL6,  IL8,  TNFAIP3,  etc.  (Tian  et  al., 
2005),  10  immune  effector  response  genes,  and  10  y-interferon  responsive  genes.  Multiple 
mediators  of  tissue  remodeling  were  also  expressed,  i.e.,  6  extracellular  matrix-degrading 
proteases,  and  4  serine  protease  inhibitors.  Conversely,  12  extracellular  matrix  constituents  were 
induced,  suggesting  that  the  cells  are  conditioned  to  reconstitute  the  basement  membrane. 
Annotation  of  the  co-culture  genes  (http://david.niaid.nih.gov/david/ease.htm)  showed  an 
enrichment  of  genes  participating  in  particular  biology.  The  upregulated  genes,  listed  in  order  of 
significance,  were  composed  of  response  to  pathogen  and  external  stimulus,  genotoxic  stress, 
inflammatory  and  innate  immune  defense,  cell-cell  interaction,  and  tissue  remodeling 
(Supplementary  Table  II).  Thirty  nine  genes  were  down-regulated  in  response  to  co-culture.  The 
down-regulated  gene  list  annotation  was  enriched  for  the  following:  cell  cross-talk,  signal 
transduction,  lipid  metabolism,  development  and  morphogenesis.  The  inflammatory  gene 
expression  changes  are  remarkable,  considering  the  absence  of  inflammatory  B-,  T-cells  or 
macrophages  in  this  simple  binary  in  vitro  culture  system.  As  4  of  the  cell  pairs  described  above 
involved  both  epithelial  and  fibroblast  cloned  cell  lines,  the  contribution  of  rare  immune  cells  in 
the  co-cultures  is  highly  unlikely. 

The  co-culture  induced  gene  signature  is  elicited  predominantly  by  soluble  factors. 

To  better  define  the  signals  that  induced  the  gene  expression  changes  in  co-culture,  we 
exposed  the  single  cell  cultures  of  each  cell  type  to  different  experimental  conditions  and 
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expression  profiled  (using  microarrays)  the  response,  aiming  to  find  conditions  that  best  mimic 
the  co-culture  response.  Gene  expression  profiles  for  each  experimental  condition  were 
compared  to  baseline  expression  of  pure  populations  of  the  same  cells  in  monoculture.  Genes 
were  selected  based  on  the  approach  described  in  Figure  IB.  The  co-culture  expression  response 
was  largely  recapitulated  when  fibroblasts  and  epithelial  cells  were  grown  on  opposite  sides  of 
transwell  plates  (Figure  3A).  These  results  imply  that  a  soluble  factor,  secreted  by  epithelial 
cells,  elicits  the  response  in  the  fibroblast  cells.  Recapitulating  the  co-culture  response  in  the 
transwell  experiment  also  means  that  the  original  observation  was  not  simply  based  on  cell  type 
ratio  bias  between  the  two  cell  mixtures  (see  supplement  1).  By  contrast,  the  other  experimental 
conditions  such  as  cell  density,  medium  exhaustion,  adherence  to  matrix,  cell  attachment  growth 
conditions,  all  generated  expression  responses  that  failed  to  mimic  the  co-culture  response. 
(Supplementary  Figure  1C,  Supplementary  Table  III).  These  results  suggested  that  the  cross  talk 
is  driven  by  soluble  but  not  contact-mediated  signals. 

Identifying  the  cell-specific  gene  expression  response  to  co-culture 

The  transwell  experiment  readily  identifies  which  gene  changes  in  the  co-culture  occur  in 
which  cell  type.  Eighty  two  of  the  125  genes  altered  in  co-culture  were  modulated  in  the  same 
direction  in  fibroblast  cells  used  in  the  transwell  experiment  (Table  1,  column  2F,  Supplementary 
Figure  1C),  whereas  58  of  the  125  genes  were  altered  in  epithelial  cells  (Table  1,  columns  2E 
and  2008).  The  differential  response  to  the  co-culture  is  also  apparent  at  the  protein  level.  While 
MMP3  is  upregulated  in  the  fibroblast,  MMP2  is  upregulated  in  the  epithelial  cell  (Figure  3B). 
The  MMP2  expression  illustrates  the  advantage  of  the  transwell  system  over  the  cell  mixture,  as 
the  fibroblast  constitutive  expression  of  the  gene  would  quench  the  induction  observed  in  the 
epithelial  cells.  Comparable  expression  change  in  the  epithelial  and  fibroblast  cell,  for  MMP2 
and  MMP3 ,  respectively,  gave  different  ratios  of  response  in  the  initial  experiment,  since  MMP2 
is  constitutively  expressed  in  the  fibroblasts,  while  MMP3  is  not.  The  contribution  of  stroma  to 
carcinoma  invasion  through  MMPs  is  well  characterized  (Jessani  et  ah,  2004).  These  results 
demonstrate  that  direct  epithelial-fibroblast  interactions,  is  sufficient  to  induce  the  stromal 
production  of  these  enzymes  in  a  manner  reminiscent  of  cancer  stroma. 

TNFa  mediates  the  co-culture  response  in  fibroblasts. 
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A  major  advantage  of  our  experimental  approach  stems  from  the  potential  to  identify  the 
basic  inductive  signal.  In  order  to  identify  the  factors  involved  in  the  cross-talk,  we  first  tested 
the  ability  of  different  cytokines  to  elicit  fibroblasts  responses  that  were  consistent  with  those 
seen  in  the  transwell  experiments  (Figure  4).  Using  ELISA  blot  analysis  of  conditioned  media 
and  comparative  expression  profiling  of  cancer  cell  cultures  between  co-culture-reactive  and 
non-reactive  lines,  we  selected  20  secreted  factors  as  candidate  mediators  of  the  epithelial  to 
fibroblast  signaling  (Supplement  2,  Supplementary  Figure  2A  and  B).  We  added  the  candidate 
cytokines  to  fibroblast  monocultures  and  followed  their  gene  expression  response  using 
quantitative  real-time  PCR  (Figure  4A,  cyan  treatments).  The  data  is  represented  in  normalized 
fold  expression  change,  as  compared  to  serum  free  media.  As  reference,  we  included  a  number 
of  the  characterized  co-culture  responses  so  that  similar  responses  would  co-cluster  (Figure  4A, 
columns  1,  2,  6-8  black  treatments).  We  found  that  the  fibroblast-epithelial  cell  co-culture 
signature  was  most  closely  recapitulated  when  the  fibroblasts  were  treated  with  epithelial 
conditioned  medium  (column  3),  or  recombinant  ILip  or  TNFa  proteins  (columns  4-5).  Not  all 
of  the  co-culture  responsive  genes  were  useful  in  identifying  the  driving  factor.  For  example, 
RAB31  and  some  of  the  adjacent  genes  are  induced  promiscuously,  while  INFy  responsive  genes 
(e.g.  STAT1,  CXCL9,  MX1,  IFITM1,  OAS1,  etc.)  are  almost  exclusively  expressed  in  response 
to  co-culture.  These  results  suggest  that  either  ILip  or  TNFa  is  sufficient  to  reproduce  the 
majority  of  the  fibroblast  response  to  the  co-culture. 

To  test  whether  ILip  or  TNFa  signaling  in  the  fibroblast  is  necessary  for  the  co-culture 
response,  we  added  neutralizing  antibodies  against  TNFa  to  the  co-culture  experiments  (as  in 
Figure  1A).  Anti-TNFa  antibodies  specifically  inhibited  the  co-culture  response,  while 
antibodies  against  TGFP  did  not  (supplementary  Figure  1C).  To  further  substantiate  this 
observation,  we  introduced  siRNA  molecules  to  the  fibroblasts,  and  tested  these  knock  down 
fibroblasts  in  co-culture  reactions  (as  in  Figure  2)  with  ovarian  cancer  cell  line  2008  (Figure  4B, 
red  treatments).  The  genes  that  were  promiscuous  were  not  included  in  this  analysis.  Knock 
down  of  essential  TNFa-signaling  components,  NIK  (IicKy,  column  17)  and  TNFR1  (column 
16),  in  fibroblasts  resulted  in  a  grossly  attenuated  co-culture  responses.  Knock-down  of  the 
IL1R1  (column  13),  MAP3K14  (column  14),  and  STAT3  (column  15)  attenuated  the  co-culture 
response  to  a  lesser  extent,  while  irrelevant  siRNA  molecules  (IL18R,  MMP26,  SMAD4  and 
CCR1)  did  not  affect  the  co-culture  response  (columns  2-5).  Since  the  ILip  gene  is  itself  co- 
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culture  responsive,  it  is  hard  to  distinguish  downstream  effects  from  a  primary  driving  role. 
However,  the  siRNA  results  with  IL1R1  suggest  that  IL 1  (3  is  indeed  a  secondary  product  of  the 
co-culture,  rather  than  the  primary  driver.  These  results  suggested  that  TNFa  or  a  TNFa-like 
molecule  secreted  by  the  epithelial  cells  drives  a  robust  inflammatory  response  in  fibroblasts. 
This  conclusion  is  further  substantiated  by  analysis  of  the  human  cancer  expression  profiles 
(supplement  4). 

Time  course  of  the  co-culture  response. 

TNFa  protein  levels  were  increased  in  co-culture  conditioned  media  compared  to  the 
mono-culture  conditioned  medium  (supplementary  figure  2A).  Since  TNFa  mRNA  is  itself  not 
co-culture  responsive  at  the  24  hour  time  point,  we  aimed  to  identify  the  time  point  at  which 
TNFa  mRNA  is  first  induced  (assuming  mRNA  accumulation  drives  the  protein  accumulation). 
We  therefore  profiled  the  time  course  dynamic  of  the  co-culture  response  for  the  2008  cell  pair 
(Figure  5A).  Genes  that  are  upregulated  or  down-regulated  at  the  24  hour  time  point  are  marked 
by  red  and  green  lines,  respectively.  The  TNFa  mRNA  level  increased  eighteen  fold  in  the  co¬ 
culture  at  the  8  hour  time  point  (Figure  5A,  bold  line).  By  24  hours,  TNFa  is  no  longer 
differentially  expressed.  TNFa  mRNA  accumulation  occurred  in  the  epithelial  cell  at  8  hours,  at 
a  time  point  preceding  the  TNFa-driven  co-culture  response  in  the  fibroblast  at  24  hours.  This  is 
consistent  with  the  assumption  of  reciprocal  dynamic  interaction  initiated  by  a  signaling  factor 
constitutively  secreted  by  the  fibroblast  cell. 

Identifying  upstream  regulators  of  TNFa  expression. 

To  investigate  the  factor(s)  responsible  for  inducing  the  the  co-culture  regulated  gene 
signature  in  the  epithelial  cells,  (especially  MCSF,  IL1  (3  and  TNFa),  we  profiled  the  response  of 
epithelial  mono-cultures  to  soluble  factors  using  real  time  PCR  (Figure  5B).  Candidate  soluble 
factors  were  based  on  genes  in  the  published  literature  that  mediate  mesenchymal  signaling  to 
the  epithelial  cell,  and  the  ELISA  results  with  fibroblast  conditioned  media  (Supplementary 
Figure  2).  TNFa  was  upregulated  228-fold  in  response  to  combined  treatment  with  either  IL8 
with  EGF  (column  3),  or  LIF  with  EGF  (column  4).  TNFa  had  induced  its  own  expression 
(column  5),  while  neither  IL8  (column6),  LIF  (column  7)  or  other  candidates  (columns  8-14) 
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could  induce  TNFa  expression  by  themselves.  The  response  of  epithelial  cells  to  TNF  involves 
an  inducible  autocrine  cascade  (Janes  et  ah,  2006).  Furthermore,  siRNA-mediated  knock  down 
of  gpl30/IL6ST  (column  15),  uPAR  (column  16)  or  SMAD4  (column  17)  did  not  reduce  TNFa 
upregulation  in  the  co-culture,  while  the  knock  down  of  STAT3  (column  18),  IFNyRl  (column 
19)  or  IL8  receptor  (column  20)  eliminated  the  TNFa  response  to  co-culture.  These  results 
demonstrate  that  IL8,  in  the  presence  of  EGF,  is  necessary  and  sufficient  to  induce  TNFa  in  the 
epithelial  cell. 

Epithelial-Fibroblast  expression  signature  is  recapitulated  in  human  cancer. 

It  is  conceivable  that  the  expression  signature  described  above,  and  the  signaling  that 
leads  to  this  coordinate  gene  expression  change,  is  unique  to  the  experimental  cell  culture  system 
in  which  it  was  defined  and  not  representative  of  the  epithelial-fibroblast  interaction  response  in 
cancer.  To  address  the  "coherence"  of  the  genes  induced  in  the  in  vitro  experiment,  in  terms  of 
their  expression  in  vivo  (i.e.  the  degree  to  which  they  are  co-expressed  in  human  cancers),  we 
examined  expression  profiles  of  174  carcinoma  profiles  (Su  et  al.,  2001).  For  this  analysis  we 
defined  a  “query”  using  the  54  genes,  (rationale  in  supplement  3)  that  were  strongly  and 
frequently  upregulation  in  co-culture  (Table  1,  bold  genes).  We  ranked  genes  according  to  how 
strongly  they  correlate  with  each  gene  of  the  query  list  (pearson  correlation),  and  calculated  the 
quartiles  and  mean  of  the  ranks  of  the  remaining  53  genes  on  the  list.  Figure  6A  is  a  density  plot 
that  draws  the  quartiles  and  mean  like  topological  curves  along  a  matrix  of  the  U95A  human 
gene  list,  where  each  position  on  the  X-axis  defines  one  of  the  54  query  genes,  according  to 
which  the  rest  of  the  U95A  genes  are  ranked  in  the  Y-axis.  The  lower  the  mean  and  quartile  are 
along  the  (Y-axis)  whole  genome  percentile,  the  more  significant  the  coexpression  of  the  overall 
54-gene  query  list  is  with  this  gene.  The  54  parallel  ranks  are  ordered  in  increasing  mean-rank 
order.  By  repeating  this  procedure  for  80,000  random  gene  lists  of  equivalent  size,  we  found  that 
the  mean  at  each  position  never  exceeded  that  of  the  co-culture  query  list  until  the  last  gene.  A 
typical  random  gene  list  distributes  the  mean  between  the  40-60%  bars  (Figure  6B).  Fifty  three 
gene  from  the  cell-cycle  module  (Iyer  et  al.,  1999)  were  even  more  tightly  co-expressed  than  the 
co-culture  genes  (Figure  6C),  only  for  the  top  60%  of  the  cell  cycle  query  list. 
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Defining  the  genes  that  constitute  the  epithelial-fibroblast  interaction  signature  in  human 
cancer. 

One  of  the  most  important  uses  of  whole-genome  expression  data  is  the  discovery  of  new 
genes  with  similar  function  to  a  given  list  of  genes  (the  query)  already  known  to  have  closely 
related  function.  Genes  that  show  tight  levels  of  co-expression  in  a  large  number  of  diverse 
microarray  experiments  are  likely  to  function  together  (Eisen  et  ah,  1998),  either  through  direct 
protein-protein  interactions  (Rhodes  et  ah,  2005),  by  directly  controlling  the  expression  of  one 
another  (Harbison  et  ah,  2004;  Long  et  ah,  2004),  or  by  participating  in  biological  processes  that 
are  coordinated  (Cam  et  ah,  2004),  possibly  occurring  in  two  neighboring  cells  types  in  the  tissue 
(Becker  et  ah,  1997).  The  binary  co-culture  in  vitro  experiment  is  an  over-simplified  model  of 
the  interactions  that  occur  in  the  tumor  microenvironment.  The  gene  recommender  algorithm 
ranks  genes  according  to  how  strongly  they  correlate  with  a  set  of  query  genes.  Using  the  same 
54  gene  query  described  above,  we  extended  our  investigations  beyond  the  binary  cross-talk  and 
the  conditions  used  in  vitro  to  tightly  linked  events  in  vivo.  We  used  the  gene  recommender 
algorithm  (Owen  et  ah,  2003),  with  the  54  genes  query  to  study  the  expression  profiles  of  three 
published  datasets  with  a  total  of  712  primary  human  tumors  (Ramaswamy  et  al.,  2001;  Su  et  al., 
2001;  Tothill  et  al.,  2005),  from  24  different  organs  and  of  different  histological  subtypes.  Gene 
recommender  output  genes,  i.e.  the  genes  that  most  strongly  correlated  with  the  in  vitro  co¬ 
culture  query  list,  are  ranked  for  each  of  the  published  datasets  independently  (Supplementary 
Table  II).  A  high  degree  of  coherence  was  observed  between  the  top  ranked  lists  in  the  in  vitro 
system  and  the  in  vivo  data.  This  response  was  also  conserved  between  independent  in  vivo 
datasets  (colored  genes).  The  gene  recommender  output  was  comparable  when  the  query  list  was 
substituted  with  genes  derived  from  two  recently  published  independent  epithelial-fibroblast  co¬ 
culture  expression  profiling  studies  (Fromigue  et  al.,  2003;  Gallagher  et  al.,  2005), (data  not 
shown).  Many  of  the  top  ranked  genes  are  bone  fide  immune  cell  markers,  such  as  CD  11c, 
CD68,  CD3,  NK4,  for  dendritic  cells,  T-lymphocytes  and  natural  killer  cells,  respectively.  This 
observation  suggests  that  the  tissue  where  the  query  genes  are  most  strongly  expressed  contained 
prevalent  immune  cell  infiltrates.  This  is  consistent  with  our  in  vitro  observation  that  the 
epithelial-fibroblast  interaction  triggers  the  expression  of  potent  chemokines  which  have  been 
shown  to  effectively  trigger  immune  cell  infiltration  (Curiel  et  al.,  2004;  Milliken  et  al.,  2002; 
Zhang  et  al.,  2003).  Also  prevalent  in  the  genes  correlated  with  the  query  are  proteases, 
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suggesting  that  the  induction  of  the  co-culture  responsive  genes  in  human  tumors  is  driven  by 
disruption  of  normal  tissue  architecture  and  basement  membrane  components.  We  show  here  that 
the  response  of  epithelial  and  fibroblast  cells  to  simple  commingling  in  in  vitro  co-culture  is 
reminiscent  of  proteolytic-enzyme  rich  tissues  in  vivo.  This  response,  in  turn,  induces  further 
proteases,  such  as  MMP1.  Based  on  these  observations,  we  propose  the  term  Tissue  Injury 
Response  (TIR)  to  the  gene  signature  we  identified  in  co-culture. 

Two  types  of  wound  healing  in  cancers. 

A  recent  attempt  to  identify  a  common  feature  in  human  cancer  expression  profiles 
suggested  that  serum  response  of  fibroblasts  recapitulates  a  consistent  aspect  of  cancer 
expression  profiles  and  further  suggested  that  this  response  also  occurs  during  wound  healing 
(Chang  et  ah,  2005;  Chang  et  ah,  2004).  However,  although  referred  to  as  a  wound  healing 
signature,  the  DAVID  annotation  for  this  published  fibroblast  response  to  serum  is  mainly 
composed  of  cell  growth  and  proliferation.  By  contrast,  the  TIR  genes  bear  no  resemblance  to 
cell  cycle  or  cell  growth  and  proliferation  signatures.  We  clustered  expression  profiles  from 
gastric  cancers  focusing  on  gene  recommender-derived  genes  of  the  two  gene  signatures;  the 
published  wound  healing  response  and  the  TIR.  Indeed,  the  two  gene  expression  signatures 
exhibited  elevated  expression  only  in  some  of  the  patient  specimens  (Figure  6C).  More  recently, 
the  wound  healing  response  was  linked  to  two  genome  aberrations  in  the  tumor  cells,  MYC  and 
CSN5  (Adler  et  al.,  2006),  which  is  in  agreement  with  the  DAVID  annotation  of  their  list,  i.e. 
representing  epithelial  cell  growth  (MYC),  rather  than  reactive  fibroblasts.  This  figure  illustrates 
that  the  genes  derived  from  the  serum  response  of  fibroblasts,  and  the  TIR  genes,  are  both  tightly 
coordinated  expression  responses,  yet  not  related  to  each  other  and  with  very  little  cross  over. 
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Discussion: 

Genomics  has  the  potential  to  revolutionize  the  diagnosis  and  management  of  cancer  by 
offering  an  unprecedented  comprehensive  view  of  the  molecular  underpinnings  of  pathology. 
Computational  analysis  is  essential  to  transform  the  masses  of  generated  data  into  a  mechanistic 
understanding  of  disease.  We  hypothesized  that  the  introduction  of  microarray  gene  expression 
measurements  from  distinct,  highly  defined  and  controlled  in  vitro  experiments  will  facilitate  the 
analysis  and  interpretation  of  human  tumor  expression  profdes,  and  will  allow  the  inference  of 
biological  mechanisms  from  expression  patterns. 

Recent  publications  produce  compelling  evidence  of  the  important  contribution  of  cancer 
microenvironment  to  cancer  pathology  and  phenotypes.  A  number  of  mechanisms  have  been 
proposed  to  explain  how  the  inflamed  stroma,  without  itself  undergoing  transformation,  may 
contribute  to  the  transformation  process  of  epithelial  cells.  When  tested  in  animal  models, 
genetic  perturbation  of  some  critical  signaling  pathways  in  the  stroma  can  lead  to  carcinoma 
(Bhowmick  et  al.,  2004;  Campisi,  2005;  Greten  et  al.,  2004;  Hagemann  et  ah,  2004;  Hiratsuka  et 
al.,  2002;  Pikarsky  et  al.,  2004;  Robinson  et  al.,  2003;  Weaver  et  al.,  2002).  The  principle  behind 
each  of  these  models  of  stromal  perturbation  is  based  on  modification  of  inflammation,  ECM 
production  or  remodeling.  For  example,  knock  out  mice  of  TNFa  and  NFicB  (essential  for 
TNFa-response)  attenuate  cancer  formation  in  animal  models  of  chronic  infection  (Greten  et  al., 
2004;  Pikarsky  et  al.,  2004).  Using  expression  profiling  of  both  experimental  models  and  human 
cancer  specimens,  we  showed  that  one  can  define  tumor-host  signaling  events  in  vitro  and 
measure  their  frequency  of  occurrence  in  cancer  patient  specimens.  In  particular,  this  work 
focuses  on  the  spontaneous  response  of  epithelial  and  fibroblast  cells  to  their  co-culture, 
irrespective  of  the  transformation  status  of  the  cells  (normal  cell  types  produce  a  similar 
response).  The  epithelial-fibroblast  interaction  in  vitro  induced  the  same  gene  groups  as  observed 
in  the  stromal  genetic  mutants  studied  in  animal  models  (inflammation,  cell  signaling  and  ECM 
production  and  remodeling).  In  addition,  siRNA-mediated  knock  down  of  critical  TNFa- 
signaling,  attenuated  the  gene  expression  provoked  by  the  co-culture.  It  has  been  suggested  that 
the  context  by  which  mutations  in  stroma  lead  to  carcinoma  is  perturbation  of  the  epithelial- 
stromal  communication  (Jacks  and  Weinberg,  2002;  Fittlepage  et  al.,  2005).  Our  results  provide 
evidence  that  supports  this  suggestion. 
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Although  inflammatory  signals  are  highly  associated  with  tissue  injury,  the  mechanisms 
of  induction  of  ILip  and  TNFa,  master  regulators  of  inflammation,  have  mainly  been  studied 
using  bacterial  wall-derived  irritants,  such  as  lipopolysaccharide.  The  high  prevalence  of 
immune  recruiting  chemokines  in  the  co-culture,  as  well  as  the  tight  association  of  the  co-culture 
signature  with  the  leukocyte  infiltrate  component  of  tumor  profiles,  suggests  that  in  vivo  immune 
cells  would  respond  to  the  chemokine  signal  by  infiltrating  into  the  corresponding  tissue.  This 
model  suggests  that  the  initiation  of  innate  immunity  and  inflammatory  response  is  governed  by 
altering  the  epithelial-fibroblast  interactions,  ostensibly  as  a  result  of  basement  membrane 
breakdown  and  wounding.  The  dynamic  reciprocal  exchange  of  factors  between  the 
mesenchymal  and  ectodennal  cell  layers  we  identified  leads  to  a  transient  increase  in  gene 
expression.  Most  of  the  identified  genes  are  no  longer  differentially  expressed  by  48  hours. 
Matrix  metalloproteinases  (MMP)  were  amongst  the  genes  most  strongly  induced  by  co-culture 
(Table,  column  1).  MMP  enzymes  have  been  implicated  in  cancer  metastasis  and  invasion 
(Lynch  and  Matrisian,  2002;  Matrisian  et  al.,  2001;  Zucker  et  al.,  2000).  We  also  found  induction 
of  active  MMP  enzymes  in  vitro  (data  not  shown).  The  proteases  probably  induce  the  transient 
response  observed  in  the  in  vitro  co-culture  by  bringing  the  epithelial  and  fibroblast  cells 
together.  Considering  the  transient  nature  of  the  response  in  vitro,  it  is  likely  that  the  elevated 
levels  of  correlated  genes  in  the  cancer  tissue  are  the  result  of  a  steady  state  of  tissue  breakdown 
and  transient  repair  response,  suggesting  that  cancer  is  like  a  chronic  wound  (Dvorak,  1986). 
Indeed  the  essence  of  cancer  invasion  is  continuous  protease  mediated  tissue  degradation.  The 
role  of  tissue  architecture  in  cancer  progression,  and  the  role  of  immune  surveillance  (especially 
in  ovarian  cancer)  are  well  established.  We  show  here  that  the  two  processes  are  coordinated. 
Consequently,  the  role  of  TNFa  in  cancer  promotion  (Luo  et  ah,  2004;  Pikarsky  et  ah,  2004), 
including  in  ovarian  cancers  (Lancaster  et  ah,  2004;  Pohl  et  ah,  2005)  is  restricted  to  chronic 
inflammation,  as  the  transient  reaction  will  not  lead  to  cancer. 

Recent  gene  knockouts  established  that  Amphiregulin  (a  ligand  of  the  EGFR  family) 
likely  provides  the  epithelial-stromal  signal  required  for  normal  mammary  ductal  morphogenesis 
(Stemlicht  et  ah,  2005;  Troyer  and  Lee,  2001)  and  protection  from  intestinal  damage  (Lee  et  ah, 
2004).  In  addition,  senescent  fibroblasts,  which  lead  to  carcinoma  (Campisi,  2005),  facilitate 
epithelial  prostate  transformation  through  Amphiregulin  (Bavik  et  ah,  2006).  Interestingly, 
senescent  fibroblasts  express  elevated  levels  of  IL8  and  MMPs  (Parrinello  et  ah,  2005), 
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suggesting  that  the  transient  response  we  observed  in  co-culture  may  be  constantly  firing  in  the 
senescent  fibroblast  irrespective  of  their  neighboring  cells.  TNFa,  IL8  and  EGFR-signaling 
pathways  all  take  part  in  the  inflammatory  response.  We  have  show  here  that  EGFR  family 
signaling  together  with  IL8  leads  to  atypical  expression  of  TNFa  by  the  epithelial  cell,  and 
consequently  induces  fibroblast  gene  expression,  typical  to  inflammation. 

Confounding  the  research  of  tumor-stroma  interaction  are  the  paradoxical  effects  that  the 
stroma  can  have  on  cancer  epithelium:  stromal  cells  are  reported  to  suppress  (Kohno  et  ah,  1982; 
Sugiyama  et  al.,  2001)  and  enhance  (Camps  et  ah,  1990;  Moore  et  ah,  1999;  Olumi  et  ah,  1999) 
malignancy.  Variations  in  human  genome  sequences  modulate  the  susceptibility  to  cancer  and 
other  diseases.  Deciphering  the  extent  of  stromal  contribution  to  cancer  progression  and  inferring 
the  relative  contribution  of  different  signaling  pathways  involved  in  communication  between  the 
stroma  and  epithelium  is  limited  by  a  lack  of  knowledge  about  the  key  locus  of  variation  among 
individual  patients’  genomes.  With  the  release  of  HapMap  (Cheung  et  al.,  2005),  we  predict  that 
haplotype  differences  in  the  constituents  that  mediate  the  co-culture  response  would  lead  to 
variation  of  the  response  through  hyperactive  or  hypoactive  signaling  (Chen  et  al.,  2006;  Guo  et 
al.,  2005;  Kang  et  al.,  2006;  Liu  et  al.,  2005;  Seifart  et  al.,  2005;  Shih  et  al.,  2006;  Zecevic  et  al., 
2006).  Therefore,  the  ability  to  correlate  the  expression  patterns  of  genes  involved  in  the  discrete 
aspects  of  tumor-stroma  interaction  to  cancer  promoting  and  attenuating  effects  may  help  to 
resolve  the  controversy.  For  example,  genes  like  CCL18  (Leung  et  al.,  2004)  and  CXCL9  (Zhang 
et  al.,  2003)  are  predictive  of  improved  patient  outcome,  whereas  CXCL12  (Oonakahara  et  al., 
2004)  and  CCL22  (Curiel  et  al.,  2004)  are  associated  with  worse  outcome.  The  co-culture 
response  of  epithelial  and  fibroblastic  cells  displayed  significant  variation  in  the  expression  of 
these  genes,  across  multiple  isolates,  in  a  manner  similar  to  the  expression  profile  variation  of 
human  cancer  specimens.  Specifically,  CXCL9  was  induced  in  the  co-cultured  epithelial  cells  in 
a  frequency  similar  to  that  observed  in  patient  specimens  by  Zhang  et  al.  (53%  in  our  assays  and 
58%  in  the  published  data).  Identifying  the  upstream  events  that  control  such  genes  may  help 
resolve  the  contention  over  the  stromal  effect  on  epithelial  cancers.  We  identified  soluble  factors 
that  were  necessary  and  sufficient  for  the  in  vitro  co-culture  response,  including  TNFa,  IL8,  and 
EGF-like  ligands,  thereby  making  these  the  first  candidates  to  test. 

The  explosion  of  information  about  the  molecular  pathways  of  cell  transformation  has  led 
to  several  highly  successful  targeted  approaches  to  cancer  treatment,  such  as  the  recent  success 
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of  imatinib  (Gleevec),  a  tyrosine  kinase  inhibitor,  for  treating  chronic  myelogenous  leukemia. 
The  “central  dogma”  of  cancer  research  is  that  mutations  in  oncogenes  and  tumor  suppressor 
genes  converge  into  a  downstream  decision  of  a  cell  whether  or  not  to  divide.  However, 
assuming  that  the  levels  of  gene  expression  indicate  the  degree  of  activity  of  a  given  biological 
program,  then  the  fact  that  many  cancer  tissue  do  not  upregulate  cell  cycle  controlled  genes 
suggests  that  cell  cycle  is  not  the  be  all  and  end  all  of  cancer.  In  particular  as  presented  in  Figure 
6,  many  gastric  cancer  specimens  do  not  express  higher  levels  of  cell  cycle  genes  than  non- 
malignant  tissues.  Thus,  not  all  cancers  depend  on  deregulated  cell  cycle,  division  and  growth. 
Rather,  some  cancers  upregulate  the  expression  of  the  genes  induced  in  the  co-culture  reaction  as 
we  described  it  (the  TIR  expression  module).  Co-regulated  expression  modules,  highly 
resembling  the  TIR  expression  module  we  observed,  have  been  described  in  a  high  throughput 
analysis  of  cancer  expression  profdes  (Segal  et  ah,  2005;  Segal  et  al.,  2004),  as  well  as  in 
ovarian,  breast  and  gastric  cancer  datasets  (Boussioutas  et  al,  in  preparation).  Recent 
publications  demonstrate  that  chemokines  (Curiel  et  al.,  2004;  Zhang  et  al.,  2003),  COL6A3 
(Sherman-Baust  et  al.,  2003;  Simojoki  et  al.,  2003),  as  well  as  the  expression  module  that 
overlaps  with  our  TIR  expression  module  (Spentzos  et  al.,  2004),  are  all  predictive  of  ovarian 
cancer  survival.  Thus,  understanding  the  factors  outside  the  cancer  cell  that  drive  it  to  invasion 
and  metastasis  or  keep  it  contained  could  provide  an  entirely  new  avenue  of  therapy.  An  avenue 
that  carcinoma  cells  may  find  harder  to  evade,  and  may  prove  relevant  to  a  larger  fraction  of 
patients  than  current  targeted  therapy  approaches,  particularly  for  patients  that  do  not  respond  to 
standard  cytotoxic  treatments  (Sherman-Baust  et  al.,  2003). 
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Experimental  Procedures 

Culturing  cells  from  patients:  For  Normal  Ovarian  Surface  Epithelial  (NOSE),  we  used  cervical  brush  scrapping  of  the  ovary,  and 
cultured  them  directly.  For  Ovarian  Stroma,  0.5x0. 5x1  cm  Ovarian  tissue  was  collected  from  consented  women,  in  DMEM  with  10%  FCS  from 
theatre  in  Monash  Medical  Centre,  and  sent  to  us  promptly,  minced  with  razor  blade,  washed  in  Trypsine-Versen,  and  then  incubated  for  2-6 
hours  in  37°C  in  10  ml  (per  0.25  gr  tissue)  of  RPMI  with  3  mg/ml  Collagenase  II  (Wortington®),  24  g/ml  DNase  I  (Boehringer  Manheim®). 
When  minced  pieces  homogenize,  the  cells  are  washed  three  times  in  PBS  with  2  %  FCS,  and  plated  in  a  culture  flask  with  fresh  medium  (Flask 
precoated  with  BD  Collagen  I®).  If  medium  contains  Gentamycin,  it  will  select  for  NOSE  cells,  if  not,  Stroma  cells  will  prevail  by  default.  This 
notion  was  confirmed  by  FACS  staining  with  anti-cytokeratin  and  Anti  Epithelial  Membrane  Antigen,  or  by  microarray  profile. 

Basic  (plastic)  Co-culture:  Stroma  cells  are  plated  in  50%  confluence  and  allowed  to  reach  80%  confluence  (24-48  hours).  Cells  are 
washed  twice  with  PBS,  and  overlayed  with  fresh  Serum  Free  Medium  (RPMI,  25mM  HEPPES,  supplemented  with  Insulin,  Transferin, 
Albumin,  and  Selenium  from  SIGMA®).  After  48  hours,  epithelial  cells  from  80%  confluence  are  trypsinized,  washed  twice  in  serum  free 
medium  plus  Soy  bean  trypsin  inhibitor,  and  plated  in  fresh  SFM,  either  on  top  of  Stroma  flask,  or  in  new  plastic  dish  (equivalent  amount  of 
Stroma  cultured  are  simply  given  fresh  SFM).  After  varying  time,  for  the  stroma  of  the  experiments  16  hours,  cells  are  trypsin-harvested,  washed, 
flash  frozen,  and  subnjected  to  RNA  preparation,  using  Triazol  (Invitrogen®)  followed  by  RNeasy  (Qiagen®).  Isolated  Stroma  and  epithelial 
cultures  are  pooled  prior  to  RNA  preparation  (E+S).  RNA  is  harvested  from  co-cultures  of  Stroma  and  epithelium  (ES). 

Conditioned  medium  experiments: 

Stromal  and  OSE  cells  were  cultured  to  >70%  confluency,  washed  twice  with  PBS,  and  growth  was  continued  for  48  hours  to  produce  the 
corresponding  conditioned  media  in  SFM.  Conditioned  media  were  filtered  through  0.2  PVDF  membrane  (WHATMAN®).  Fresh  cells, 
Stromal  or  epithelial  were  plated  and  serum  starved  as  above,  and  then  fresh  SFM  or  SFM  supplemented  with  5  fold  dilution  of  their  own  or  their 
counterpart  cell  CM,  or  with  commercial  recombinant  TGF  ,  bFGF  (Sigma  ®)  or  TNF  (BD®)  in  20  ng/ml  concentration.  Cells  grow  in  the 
new  medium  for  24  hours,  harvested  and  processed  as  above. 

Contact  dependent  cues  assays: 

To  identify  co-culture  regulated  genes  that  are  direct  membrane  protein  interactions-dependent,  we  performed  the  co-culture  experiment,  using 
tissue  culture  inserts  (0.4  M  pores  NC,  Coming®)  with  Stromal  cells  first  plated  on  upper  side  of  the  insert,  and  then  Epithelial  cells  plated 
either  on  the  same  side  (contact  allowed),  or  on  the  bottom  side  (only  soluble  factors  allowed).  The  two  co-cultures  are  directly  compared  to  each 
other  in  a  microarray  experiment. 

ECM  signal  assays: 

Collagen  I,  Matrigel,  human  normal  ECM  (all  from  BD®)  were  added  to  SFM  in  50  g/ml  concentration  and  used  to  overlay  T175  culture  flasks 
(Coming®  filter  cap)  overnight.  The  next  day,  the  flasks  were  washed  with  PBS  and  OSE  cells  were  plated  and  grown  in  either  fresh  un-coated 
flasks,  or  the  coated  flasks.  After  24  of  growth,  cells  were  harvested  and  processed  as  above. 

Organotypic  cultures: 

Rat  tail  Collagen  I  (BD®)  is  mixed  with  50%  serum  50%  RPMI  medium,  neutralized  and  allowed  to  gel  on  the  bottom  of  a  tissue  culture  insert 
(0.4p  Fibrobalsts  are  mixed  with  Rat  tail  Collagen  I,  neutralized  and  plated  on  culture  dish  inserts.  Stromal  cells  plated  in  collagen  1  (Carter  et  al., 
1990)  proliferate  and  metabolise  the  extracellular  matrix.  We  allowed  the  Stroma  to  grow  for  a  week,  while  the  gel  shrinks.  Epithelial 
Arrays:  cDNA  arrays  were  processed  as  published(Tothill  et  al.,  2005).  Affymetrix®  U133+2  arrays  were  processed  according  to  manufacturers’ 
instructions. 

Real  time  PCR: 

Twenty  seven  genes  were  scored  by  both  primer  sets  designed  in  house,  using  SYBR-green  as  well  as  TaqMan™  primer  sets  from  ABI®.  The 
independent  primer  sets  showed  average  correlation  of  R=0.85  with  each  other  and  R=0.71  with  the  array  results. 

The  response  to  cell  density  was  measured  by  comparing  cultures  of  25%  and  90%  confluence.  The  response  to  medium  Exhaustion  was  tested 
by  comparing  cultures  grown  in  media  from  48  hour  cultures,  with  fresh  medium.  The  response  to  adherence  to  extra  cellular  matrix  components 
was  tested  by  comparing  cells  grown  on  plastic  flasks,  with  cells  grown  on  top  of  flasks,  which  were  pre-coated  with  collagen,  matrigel,  and  fresh 
extracellular  matrix  preparations  from  the  fibroblasts  we  used  for  co-cultures.  The  response  to  attachment  altogether  was  tested  by  growing  cells 
in  RGD  peptides  with  0.1  mM  EDTA  in  medium,  with  cells  plated  on  the  same  time  into  normal  growth  conditions,  and  by  plating  cell  onto  agar 
coated  plates. 
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Supplement  1:  Methods  of  confirming  that  mixed  cell  population  compared  on  expression 
profiling  are  not  systematically  biased  towards  one  cell  type. 

Specific  detection  of  co-culture  regulated  genes  in  the  experiment  depicted  in  Figure  IB 
relies  on  there  being  equivalent  numbers  of  each  cell  type  in  the  co-culture  compared  with  the 
pooled  monoculture.  Due  to  the  robust  expression  differences  between  epithelial  and  fibroblast 
cells,  small  differences  in  epithelial  or  fibroblast  cell  number  between  the  co-culture  and  the 
pooled  monoculture  would  be  expected  to  result  in  lineage  specific  genes  being  incorrectly 
identified  as  co-culture  regulated.  To  assess  whether  cell  bias  contributed  to  our  findings,  we 
performed  FACS  analysis  with  epithelial-specific  marker,  CK7.  This  demonstrated  that 
equivalent  numbers  of  each  cell  type  were  present  in  the  samples  used  for  RNA  preparation 
(Supplementary  Figure  1A).  Furthermore,  we  reasoned  that  if  the  genes  which  are  apparently 
induced  by  co-culture  were  an  artifact  of  cell  bias,  then  all  co-culture  genes  should  simply  be 
epithelial  or  fibroblast  specific  (depending  upon  which  cell  type  predominated).  We  therefore 
competitively  hybridized  RNA  derived  from  epithelial  and  fibroblast  cultures  to  identify  cell 
specific  gene  expression  (Supplementary  Figure  IB)  and  plotted  log  ratio  values  against  average 
log  intensity  (M  versus  A  plot;  (Dudoit  et  al.,  2002)).  Genes  that  were  most  differentially 
expressed  in  co-culture  were  highlighted  (Supplementary  Figure  IB,  red  triangles).  Although 
many  of  the  genes  induced  were  normally  more  strongly  expressed  in  fibroblasts  they  were 
surrounded  by  a  cloud  of  fibroblast  specific  genes  that  do  not  change.  Conversely,  a  number  of 
induced  genes  are  more  strongly  expressed  in  epithelial  cells.  These  conclusions  are  supported 
by  the  recapitulation  of  co-culture  response  in  the  transwell  experiments,  which  compare 
homogenous  cell  populations. 

Supplement  2:  Choosing  candidate  factors  for  Figures  5  and  6. 

We  assumed  the  factor  that  induces  this  response  in  the  fibroblast  would  be  expressed  by  the 
epithelial  cell.  Candidate  factors  were  first  chosen  based  on  differential  expression  among  the 
two  cell  types.  TNFa  mRNA  was  four  fold  more  epithelial  specific,  so  was  TGFP3,  and 
InhibinA.  To  identify  additional  candidates,  we  harvested  conditioned  media  from  epithelial  and 
fibroblast  cells  and  examined  proteins  that  accumulate  in  those  media  using  ELISA  assays  on 
Chemicon®  human  cytokine  antibody  arrays  (Chemiarrays™  V,  VI  and  VII)  against  120 
cytokines,  and  growth  factors.  We  detected  IL6,  IL8,  MCP1  in  the  fibroblast  cells,  and  EGF,  and 
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TNFa  in  the  2008  ovarian  epithelial  cancer  cells  (Supplementary  Figure  2A).  An  epithelial 
secreted  molecule  can  also  enzymatically  activate  a  cytokine  that  is  expressed  by  the  fibroblasts. 
For  example,  MMP26  and  MMP7,  both  known  to  activate  TNFa  are  64  and  16  fold  up-regulated 
by  the  epithelial  cell.  Typical  TNF  response  genes  were  induced  in  the  co-culture  and  the 
transwell  experiments,  but  not  in  the  epithelial  conditioned  media  treatment  of  fibroblasts. 
Conversely,  the  genes  elicited  by  the  transwell  in  epithelial  cells  are  typical  IFNy  response  genes. 
This  suggests  a  dynamic  reciprocal  interaction  occurs  in  the  co-culture  and  transwell 
experiments.  We  therefore  analysed  the  factors  accumulating  in  mono-cultures  as  well  as  co¬ 
cultures  (Supplementary  Figure  2A).  Consistent  with  the  expression  array,  OPG/TNFRSF11  and 
FGF4  were  reduced  in  co-culture,  while  TNFa  protein  was  only  detected  in  co-culture  media. 
Lastly,  to  identify  additional  candidates,  we  screened  for  genes  that  are  specifically  expressed  in 
cell  lines  that  evoke  a  more  robust  co-culture  reaction.  We  used  HG_U133jplus_2  Affymetrix 
arrays  to  expression  profile  compare  A2780,  SKOV3,  DOV13,  OVCAR3  (none  responders), 
with  2008  and  OVHS  (Supplementary  Figure  2B).  For  example,  1118  and  ICE/Caspl  (activates 
ILip)  were  700  fold  differentially  expressed  among  cell  lines  2008,  OVHS,  which  evoke  the  co¬ 
culture  response,  and  cell  lines  OVCAR3,  DOV13,  and  A2780  which  do  not  respond  to  the  co¬ 
culture.  Responders  expressed  significantly  higher  levels  of  Amphiregulin,  Epiregulin  (EGF-like 
ErbB  ligand),  MAP4K5,  NIK,  and  MLK3  (signaling  kinases  activating  NFkB),  than  non¬ 
responder  cell  lines.  By  contrast,  a  co-culture  induced  gene  that  may  serve  to  insulate  the 
epithelial-mesenchymal  interaction  in  vivo  (in  the  basement  membrane)  CSPG2  was  expressed 
40-80  fold  lower  in  the  co-culture  reactive  cell  lines.  This  could  also  reflect  the  activity  of  an 
autocrine  signaling,  such  as  WNT-pathway.  Additional  criterion  to  select  candidate  factors  was 
provided  by  the  in  vivo  cancer  specimen  expression  profiles.  Secreted  factor  genes  that  are 
consistently  correlated  with  the  co-culture  response  in  the  patient  profiles  (Supplementary 
Tablell)  were  considered  candidates  and  tested  (BMP1,  INHBA,  CCL2,  MCSF,  CCL5  and 
uPA). 

Supplement  3:  Rational  behind  the  choice  of  54  genes  for  the  gene  recommender  algorithm. 

The  expression  response  to  co-culture,  as  most  expression  responses,  included  genes  that 
were  upregulated  as  well  as  genes  that  were  down-regulated.  In  seeking  genes  that  are  consistent 
with  the  overall  gene  signature,  it  is  complicated  to  include  both  upregulated  as  well  as  the 
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down-regulated  genes,  as  these  would  be  inversely  correlated  with  each  other.  Therefore,  we 
focused  our  analysis  on  “bait”  gene  list  that  was  composed  of  only  the  genes  that  were 
upregulated,  and  robustly  upregulated.  When  comparing  our  work  with  cDNA  arrays,  it  became 
apparent  that  some  genes,  such  as  MMP3  are  robust  co-culture  responders,  yet,  in  the  human 
cancer  expression  data,  these  genes  are  virtually  only  expressed  in  the  tumors,  which  complicates 
the  numerical  analysis  of  Cy3  to  Cy5  ratios.  Such  genes  were  also  excluded. 

Supplement  4:  In  vivo  expression  patterns  of  the  co-culture  responsive  genes  are  correlated 
with  TNFa  response  genes. 

Recently,  the  physiological  relevance  of  a  transcriptional  signatures,  defined  initially  in 
cells  in  vitro,  was  tested  by  analysis  microarray  data  from  hundreds  of  tumor  specimens  (Lamb 
et  ah,  2003).  For  this  analysis  we  used  the  co-culture  responsive  gene  list  described  in 
supplement  3 .  A  pattern  of  differential  expression  of  any  given  gene  across  multiple  carcinoma 
profiles,  is  used  for  ranking  genes  according  to  their  corresponding  expression  correlation  to  one 
query  gene,  such  as  the  one  shown  in  Supplementary  Figure  2C.  To  substantiate  our  conclusion 
that  TNFa-like  signal  is  responsible  for  the  majority  of  the  co-culture  reaction,  we  analyzed  the 
in  vivo  cancer  specimen  profiles  as  in  Figure  3A  and  3B  and  Supplemental  Table  II 
(Supplementary  Figure  2D).  To  support  our  conclusion  that  a  TNFa-like  signal  is  responsible  for 
the  majority  of  the  co-culture  reaction,  we  analyzed  the  in  vivo  cancer  specimen  profiles  in  the 
same  way  as  we  did  in  the  analysis  leading  to  Figures  3A  and  3B  and  Supplemental  Table  II 
(Supplementary  Figure  2D).  Our  hypothesis  now  is  that  co-culture  genes  as  a  whole  are  much 
more  strongly  associated  with  TNFa-response  genes  than  TGFP-response  genes.  Proceeding  as 
we  did  in  our  discussion  above  of  the  co-culture  gene  set  in  relation  to  the  226  ovarian  cancer 
specimens,  here  we  used  the  229  tumor  set  (Tothill  et  al.,  2005)  and  four  rankings  of  genes  on 
the  array  according  to  their  correlation  with  2  TNFa-responsive  genes,  TNFAIP3  and  DUSP5, 
and  with  2  TGFP-responsive  genes,  TGFBI  and  TGFB1I4  (Supplementary  Figures  C  and  D). 
The  selection  of  these  representative  genes  was  based  on  our  own  expression  profiling  of  TNFa 
and  TGFP  responses  (Supp  Fig  1C  and  Supp  Table  III)  using  anova  between  the  two  profiles. 
We  then  commuted  the  Wilcoxon  rank  sum  statistic  comparing  the  rankings  of  the  co-culture 
gene  sets  with  the  non  co-culture  genes,  obtaining  p-values  of  10A-13,  10A-6,  0.007  and  0.002  for 
the  above  respectively,  for  the  above  4  genes.  These  results  support  our  hypothesis  that  the  co- 
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culture  genes  as  a  whole  are  much  more  strongly  associated  with  TNFa-response  genes  than 
TGFP-response  genes.  We  scored  for  complex  correlations  across  sixty  three  independent 
treatments,  and  listed  the  treatments  that  are  most  correlated  with  each  co-culture  cell  pairs 
(Figure  2  color  of  treeview  branches,  and  see  Supplemental  Tablel  on  top  for  correlation 
coefficients).  Twelve  cell  pairs  were  significantly  correlated  with  treatments  of  fibroblast 
cultures  with  TNFa  (red),  or  MCSF1  (yellow),  while  only  three  correlated  with  EGF+IL8 
treatments  of  the  epithelial  cells  (blue).  ILip  scored  highly  close  to  TNFa  and  more  surprisingly, 
TGFP  scored  close  to  the  MCSF. 
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Figure  Legends 


Figure  1:  A.  Co-culture  versus  mono-cultures:  Equal  numbers  of  cells  were  grown  either 
in  co-culture  or  in  mono-cultures  for  16-24  hours.  These  were  harvested,  epithelial  and 
fibroblast  mono-cultures  pooled.  RNA  from  both  mixes  is  compared  either  using 
microarrays,  or  real  time  PCR.  B.  Volcano  plot  (LOD  score)  of  co-culture  response 
across  seven  repeats  of  the  experiment.  X-axis  is  M  value,  and  y-axis  is  the  B  statistic  of 
this  expression  change.  Each  dot  represents  an  individual  gene,  while  dot  color  marks  the 
different  groups  of  genes,  as  divided  by  the  statistical  properties  of  the  response  to  co¬ 
culture  in  seven  experimental  repeats.  B  statistics  represent  an  empirical  Bayesian 
statistic  (Lonnstedt  and  Speed,  2003),  T  the  student  statistic,  M  the  -log  (ratio)),  and  A 
the  log  geometric  mean(Dudoit  et  al.,  2002).  Using  these  statistical  measures  of 
significance,  we  chose  genes  that  are  part  of  the  response,  even  though  their  expression 
fold  change  was  less  than  the  usual  2  fold.  Included  genes  were  double  positive  for  either 
two  of  the  M,  B,  and  T  statistics  (empirical  threshold  as  dashed  line). 

Figure  2:  A.  Co-culture  versus  isolated  cells  for  15  independent  Ovarian  Cancer  cell 
lines.  Experiment  described  in  Figure  1A  was  repeated  with  fifteen  ovarian  cell  lines 
listed  in  column  bottom,  and  expression  profiling  was  conducted  using  real  time  PCR 
with  125  primer  sets  (Supplementary  Table  I).  Shown  are  the  co-culture  over  the  mixed 
monoculture  fold  change  for  each  cell  pair,  for  75  changing  genes,  and  nine  control 
genes.  The  data  was  normalized  relative  to  three  control  genes,  18S  rRNA,  CystA,  and 
GAPDH,  chosen  and  as  in  (Vandesompele  et  al.,  2002).  The  tree  view  of  the  different 
expression  profiles  as  organized  by  hierarchical  clustering.  Color  of  tree  view  sample 
branch  indicates  the  fibroblast  treatment-response  that  most  correlated  with  each  co¬ 
culture  cell  pair  response  (see  supplementary  Table  I,  figures  4  and  5);  red  is  TNFa/IL(3 
(R=0.72),  blue  is  InhibinA  (R=0.4),  yellow  is  TGFp/MCSF  (R=0.68).  These  correlations 
were  identified  via  a  feature  of  Genespring  GX®.  B.  Control  genes,  analyzed  as  in  A.  C. 
Box  plot  representation  of  the  overall  response  to  co-culture  in  a  logarithmic  scale.  Boxes 
represent  90%  of  the  genes.  A  bar  across  the  samples  is  on  the  level  of  two  fold  response. 
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Figure  3:  A.  An  alternative  approach  to  co-culture  that  enables  monitoring  the  response 
to  co-culture  of  each  cell  type  individually.  Epithelial  and  fibroblast  cells  from  opposite 
sides  of  a  transwell  are  compared  to  the  same  cells  growing  alone.  B.  Anti-MMP 
immunoblot  characterization  of  the  co-culture  response.  Cells  grown  as  in  Figure  3A 
(transwell  experiment)  were  collected  into  Trizol  ®  for  the  microarray  analysis.  Protein 
fraction  of  the  epithelial  and  fibroblast  cell  cultures  were  dialyzed  according  to 
manufacturer’s  instructions,  and  applied  on  SDS-acrylamide  gel  immunoblots. 

Figure  4:  A.  Identification  of  signaling  molecules  responsible  for  the  co-culture  response 
in  the  fibroblasts,  using  RT-PCR  based  profiling.  Hierarchical  clustering  compared  three 
type  of  experiments;  co-cultures  (black  labels,  as  in  Figure  2,  columns  1,  2,  6-8), 
fibroblast  response  to  epithelial  factors  (column  3),  and  fibroblast  response  to  candidate 
commercial  proteins  (cyan  labels,  columns  4,  5,  9-25).  B.  The  effect  of  single  gene 
knock-down  on  the  co-culture  response.  Fibroblasts  were  treated  with  siRNA  gene- 
specific  reagents  (red  labels),  and  then  co-cultured  with  2008  cells,  as  in  Figure  2. 
Hierarchical  clustering  compared  five  types  of  experiments;  co-cultures  (columns  1,6,9, 
10,  12),  fibroblast  response  to  epithelial  factors  (column  11),  fibroblast  response  to 
rhTNFa  commercial  protein  (column  7)  and  fibroblast  response  to  co-culture,  following 
commercial  synthetic  siRNA-mediated  knock  down  of  control  genes  (columns  2-5),  and 
genes  that  are  critical  for  the  co-culture  signaling  (columns  13-17). 

Figure  5:  A.  Time  course  experiment  of  the  co-culture  (Iyer  et  ah,  1999),  with  co-culture 
substituting  the  serum  treatment,  as  in  Figure  1A.  X-axis  describes  culture  duration  until 
harvest.  Gene  expression  lines  are  colored  according  to  their  corresponding  kinetic 
feature.  Gene  depicted  in  bold  black  is  TNFa.  B.  Identification  of  signaling  molecules 
responsible  for  the  co-culture-induced  expression  of  TNFa  niRNA  in  the  2008  epithelial 
cells,  using  RT-PCR  based  profiling.  Hierarchical  clustering  compared  four  type  of 
experiments;  co-cultures  (as  in  Figure  3,  columns  1  and  14),  co-cultures  of  siRNA 
modified  epithelial  cells  (columns  2-4,  15  and  19),  epithelial  response  to  fibroblast 
factors  (column  5),  and  epithelial  response  to  candidate  commercial  proteins  (columns  6- 
13,  16-18). 
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Figure  6:  A.  A  matrix  of  Co-expression  of  the  co-culture  responsive.  All  genes  were 
ordered  (along  the  Y-axis)  according  to  expression  correlation  with  each  of  the  54  genes 
(along  the  X-axis).  The  mean  (blue  line)  and  quartile  (top  25%  in  dashed  green,  50%  red, 
and  lowest  25%  in  dashed  purple)  rank  value  distributions  of  54  typical  co-culture  genes 
are  presented.  B.  Same  analysis  on  a  typical  random  54  gene  list.  Results  vary  for  80,000 
random  lists,  but  never  exceed  in  mean  correlation  neither  the  co-culture,  nor  the  cell 
cycle  gene  lists.  C.  The  same  analysis  on  53  cell  cycle  genes,  invoked  in  fibroblasts  upon 
serum  simulation.  D.  Hierarchical  clustering  of  43  premalignant  and  64  cancerous  Gastric 
specimens  on  260  genes,  130  representing  the  co-culture  correlated  genes  (from 
Supplemental  table  II),  and  130  genes  selected  by  serum  treatment  of  fibroblasts. 


Supplementary  Figure  1:  A.  Flow  cytometry  of  cell  mixtures  of  epithelial  and  fibroblast 
cells,  either  following  co-cultures  (left),  of  mixed  mono-cultures  upon  harvest  (right). 
The  epithelial  fraction  in  those  two  mixtures  was  37±2%  and  39±2%,  correspondingly, 
across  five  repeats.  Such  variation  would  not  justify  a  cell  type  specific  gene  showing  an 
expression  difference  higher  than  two  fold  (based  on  the  differential  expression  described 
in  F.).  B.  Mono-cultures  of  epithelial  and  fibroblast  cells  were  compared  by  expression 
profiling  (five  repeats).  Each  gene  is  assigned  an  average  M  and  A  as  above.  Top  and 
bottom  of  the  graph  represent  epithelial  and  fibroblast  specific  genes,  correspondingly. 
Red  triangles  mark  the  co-culture  regulated  genes  listed  in  Table  1.  C.  Hierarchical 
clustering  comparing  22  distinct  treatments  in  attempt  to  identify  the  cellular  event  that 
best  mimics  the  co-culture,  using  1400  genes  selected  from  the  10.5  K  in  house  cDNA 
array(Boussioutas  et  ah,  2003;  Tothill  et  ah,  2005).  These  genes  are  the  total  of  genes 
changing  in  any  of  the  experimental  treatments. 

Supplementary  Figure  2:  A.  Factors  secreted  to  the  medium  of  each  cell  line,  as 
indicated  in  the  top  right  comers,  were  identified  using  chemiarray™  ELISA  blots.  B. 
Hierarchical  clustering  comparison  of  co-culture  reactive  ovarian  cell  lines  (2008  and 
OVHS),  with  cell  lines  that  are  relatively  non-reactive  in  co-culture.  189  genes  were 
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selected  from  the  U133+2  full  genome  array  using  ANOVA  analysis.  C.  Expression 
pattern  of  a  single  gene,  TNFAIP3,  across  282  carcinoma  expression  profiles(Tothill  et 
al„  2005).  D.  Correlation  coefficient  plot  of  all  genes  to  the  expression  pattern  of  single 
gene  across  282  carcinoma  profiles  (as  in  A).  X-axis  Rank  order  of  all  genes  on  the  array. 
Y-axis  describes  the  Mann  Whitney  score,  based  on  expression  correlation  to  the  single 
gene  labeled  on  the  top  left  (blue  S-curve).  Co-culture  responsive  genes  are  marked  as 
red  circles.  P-values  are  listed  on  the  top  right  of  each  graph. 
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